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NOMENCLATURE 

A  Parameter,  A  ■  E  (1  -  T  /T 

c  s 

a  Coefficient  in  burning  rate  law,  m  *  apn 
B  Parameter,  Eq.  70,  82,  or  108 
b  Coefficient  in  Eq.  62 
c ,0^  Specific  heat  of  solid,  gas 

E  E  «  E  /R  T 
s  os 

Eg  Activation  energy  for  surface  reaction 
H  K  -  Q  /c  (T  -T  ) 

8  SC 

H  See  Eq.  85 

P 

HgT  See  Eq.  109 
h  Enthalpy 

h,  ,  h  Reference  values  of  enthalpy 
i  r  8* 

K  Eq.  128 

kc>  kg  Thermal  conductivity  of  solid,  gas 
kA  kA  "  l'/L 


iv 
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Lg  Lg  ■  hg+  -  hg_,  positive  for  endothermic  reaction 

i  t  ■  k  /me 
g  P 

Mach  number  of  product  gases 

m  Mass  flux  (gm/cm2  -  sec) 
a  Index  in  m  ■  apn 
ng  Index  in  Eq.  16 

p  Pressure 

Heat  release  per  unit  mass  in  gas  phase 

Q  Q  ~  L 
8  s  s 

q  Heat  transfer  (ergs/cm-sec) 

qf'  Heat  transfer  fluctuation  at  instantaneous  position  of 
flame  (outer  edge) 

qr'  Heat  transfer  fluctuation  at  mean  position  of  flame  (cuter 
edge) 

R  Response  function,  R  »  (m’/m)/(p  Vp) 

Rq  Universal  gas  constant 

r  Linear  burning  rate  (cm/sec) 
s  Index  in  Sq.  64,  s  -  A  in  linear  calculations 
T  Temperature 

Tc  Temperature  of  cold  solid  (x  •+  -  «°) 

t  Dimensionless  time  in  Section  5;  scale  1b  p  k  /m2e 

s  g  p 

Tf’  Temperature  fluctuation  at  instantaneous  position  of  flame 
(outer  edge) 

Tj-’  Temperature  fluctuation  at  mean  position  of  flame  (outer  edge) 


v 


V  Coefficient  in  Zq.  100 
VT  Eq.  104 

W  Eq.  65 

W0  Value  of  U  for  $  i  0,  Eq.  81 
WgT  Eq.  103 
X  Eq.  105 

Fluctuation  of  the  outer  boundary  of  the  flame 

x.,  x-  Inner  (upstream)  and  outer  (downstream)  boundaries  of  the  flame 
zone 

xg  Position  of  solid-gas  Interface 
8  Eq.  80 


a  Index  in  Eq.  16 
s 

qj  (o)  (o) 

“l  “l  “  V<0>  "  +  A2  r—  a  +  Tj  -  (2  -  5f)  A2u ' (o> 

me  T  w 

P  8 

see  Eq.  133 

o 

8  Exponent  in  Cj,  *  5^  »  Eq.  73 


* 
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,*(o) 


P„.  P„ 

c  g 


6^  *  -A2  - A2m*  0  ,  see  Eq.  1!  ’ 

w 

e  Smell  parameter,  Eq.  122 
e  Reaction  rate,  Eq.  19 
C  C  ■  exp  [J  (mcp/kg)dxl 

6  6  *  (xf  -  %)U 

Thermal  dlffusivity  of  solid,  <c  *  kc/pcc 

k  Thermal  dlffusivity  of  gae,  ic  *  k  /p  c 
8  8*  S  8®  P 

A2  Eigenvalue,  Eq.  34 

A  Complex  function  of  ft,  Eq.  4 

Ar,  A^  Real  and  imaginary  parts  of  A 

A+,  A_  Defined  In  Eq.  125 

y  Normalized  mass  flux,  v  *  m/m 


C  € 


■  I  (“A5 


dx 


5  -  (m/k)x  «  (r/<  )x 

c  c  c  c 


Density  of  solid,  gas 
t  Normalized  temperature,  t  »  T/T 


Average  temperature  introduced  in  Eq.  134 

—2 

0  Normalized  frequency,  fl  -  <cm/r 


u>  Real  angular  frequency 
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SUBSCRIPTS 


o+  Evaluated  on  the  solid  (-)  or  gas  (+)  side  of  the  origin 
s+  Evaluated  on  the  solid  (-)  or  gas  (+)  side  of  the  interface 
i  Imaginary  part 
r  R'-al  part 

f  Evaluated  at  the  instantaneous  position  of  the  flame 
f  Evaluated  at  the  mean  position  of  the  flame 

SUPERSCRIPTS 

Mean  value 
'  Fluctuation  value 

(o),  (1)  Zeroth  and  first  order  terms  in  the  expansions  (Eq.  132) 
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I.  INTRODUCTION 


It  appears  at  the  present  time  that  much  of  the  unsteady  behavior 
in  solid-propellant  rocket  chambers  can  be  traced  to  the  interaction 
between  pressure  fluctuations  and  the  burning  in  ..he  neighborhood  of  the 
solid  surface.  Thus,  an  understanding  of  the  processes  in  this  region 
when  the  pressure  is  varying  in  time  is  essential.  Ultimately,  experi¬ 
mental  results  are  required,  but  their  interpretation  rests  on  analysis, 
which  is  necessarily  approximate. 

The  various  analyses  which  have  appeared  differ  substantially  only 
in  the  model  adopted  for  the  gas  phase.  But  even  in  that  respect,  the 
basic  assumptions  involved  are  so  strong  that  almost  all  the  final  re¬ 
sults  for  the  response  to  harmonic  pressure  oscillations  have  the  same 
form,  a  two-parameter  function  of  dimensionless  frequency  (Ref.  1), 
shown  as  Eq.  71  here.  In  fact,  only  one  of  those  parameters  (denoted  by 
B)  depends  on  the  model  chosen  for  the  gas  phase;  the  other,  .  - ,  being 
associated  with  the  thermal  wave  in  the  solid  and  the  surface  reaction 
for  transformation  of  solid  to  gas.  This  conclusion  is  a  consequence 
of  four  basic  assumptions:  (1)  the  problem  is  linear  and  treated  in 
one-d in en6 iona 1  form;  (2)  the  solid  is  homogeneous  and  nonreacting,  with 
constant  properties;  (3)  the  surface  reaction  is  independent  of  pressure; 
and  (4)  the  gas  phase  responds  quasi-statically  to  changes  of  pressure. 

It  is  easy  to  relax  assumption  (3) ,  in  which  case  the  form  of  the  re¬ 
sponse  function  is  changed  only  slightly  and  one  additional  parameter  is 

introduced,  the  index  n  in  Eq.  71. 

s 

Differences  among  the  explicit  formulas  for  the  parameter  B,  ex¬ 
pressed  in  terms  of  the  properties  of  the  gas  phase,  arise  from  differ¬ 
ences  in  the  models  chosen  for  the  structure  of  the  flame  in  the  gas 
phase.  Hence,  the  detailed  interpretation  of  numerical  results  and 
experimental  data  must  differ  among  the  calculations.  So  long  as  the 
assumption  of  quasi-static  behavior  is  correctly  interpreted  and  U6ed, 
disagreement  of  this  sort  is  in  fact  associated  with  dissimilarities 
among  the  models  used  for  steady  burning.  The  major  portion  of  the 
discussion  here  is  grounded  on  one  simple,  yet  fairly  realistic,  model, 
based  on  uniformly  distributed  combustion,  which  has  been  solved  approx¬ 
imately  and  used  previously  in  a  calculation  of  the  response  function 
(Ref.  3).  The  approximations  used  in  Ref.  2,  in  addition  to  those  noted 
above,  imply  that  the  fluctuations  of  total  energy  release  in  the  gas 
phase  are  always  in  phase  with  the  pressure  fluctuations,  a  result  which 
is  open  to  question.  Alternatives  are  possible,  leading,  of  course,  to 
different  expressions  for  the  parameter  B.  Although  part  of  the 
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approximate  nature  of  the  analysis  of  Ref.  2  can  be  eliminated,  as  al¬ 
ready  noted  in  Ref.  3  and  explicitly  shown  here,  it  is  not  possible  to 
obtain  a  complete  and  unique  solution  for  the  quasi-static  problem. 

There  are  two  purposes  of  the  work  teported  here.  (1)  within  the 
feu”  crucial  assumptions  noted  above,  several  of  the  additional  approxi¬ 
mations  used  in  Ref.  2  are  examined,  and  means  of  correcting  them  are 
offered.  (2)  Two  of  the  four  basic  assumptions  listed  are  relaxed  some¬ 
what  and  the  corresponding  response  functions  are  computed:  a  separate 
calculation  is  given  for  nonquasi-static  behavior  which  does  not.  contain 
a  pressure  dependence  in  the  surface  reaction.  In  view  of  the  length  of 
this  report,  it  is  perhaps  helpful  to  summarize  the  main  points  to  be  made. 

As  in  Ref.  1,  the  calculation  of  the  response  function  is  split  into 
three  parts:  analysis  of  the  solid  phase,  construction  cf  important  con¬ 
servation  relations  for  the  sclid-gas  interface,  and  analysis  of  the  gas 
phase.  The  conservation  of  energy  applied  tc  the  interface  provides  an 
important  equation  which  is  used  as  the.  primary  relation  into  which  the 
other  results  are  8ubatltv~»d.  In  particular,  it  is  the  heat  transfer 
at  the  surface  which  must  be  found  from  the  calculations  for  the  solid 
and  gat>  phases. 

The  simplest  model  for  the  solid  will  be  used:  homogeneous,  ncn- 
reacting  material  having  constant  properties  (see  Ref.  4  and  5  for  analy¬ 
ses  in  which  chemical  activity  in  the  solid  phase  is  examined) .  Thus , 
the  treatment  of  both  the  solid  phase  and  the  interfacial  region  is 
essentially  the  same  as  in  most  other  researches  on  this  problem.  The 
necessary  information  is  collected  in  Section  Tl. 

By  far  the  most  difficult  part  of  the  problem  is  analysis  of  the 
gas  phase;  even  the  steady-state  problem  remains  at  the  present  time 
essentially  unsolved.  As  a  basis  for  all  the  work  covered  here,  the 
model  of  uniformly  distributed  combustion  is  used.  This  can  be  solved 
"exactly,"  although  certain  minor  assumptions  must  be  made.  For  example, 
if  the  combustion  zone  extends  from  some  plane  removed  from  the  inter¬ 
face  to  the  outer  boundary,  an  additional  condition,  such  as  an  "ignition 
temperature"  or  a  minimum  reaction  rate,  must  be  used  to  define  the  inner 
boundary.  However,  the  solution  can  be  found,  culminating  in  a  trans¬ 
cendental  equation  (Eq.  49)  for  the  flame  thickness.  The  formal  solution 
to  the  steady-state  problem  is  covered  in  Section  III.  The  linear  burn¬ 
ing  rate  caa  be  found  according  to  the  procedure  outlined  there,  although 
detailed  results  are  not  Included  here.* 

Perhaps  the  most  Important  assumption  of  those  noted  above  is  that 
the  gas  phase  responds  in  a  quasi-static  manner.  This  effects  enormous 


*  "Linear  Burning  Rates  for  Solid  Propellants  Having  Unifannly  Dis¬ 
tributed  Combustion,"  by  F.  E.  C.  Culick  and  G.  H.  Dehor ity  (in  prepara¬ 
tion)  . 
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simplifications  in  the  analysis.  With  this  assumption,  the  necessary 
results  can  be  obtained  in  two  ways:  either  directly ,  from  a  steady- 
state  solution,  as  discussed  in  Section  ZV-A;  or  by  solving  the  rele¬ 
vant  linearized  differential  equations,  as  covered  in  Section  IV-D.  The 
first  approach  is  simpler,  but  the  second  must  be  used  in  calculations 
not  based  on  the  assumption  of  quasi-static  behavior,  a  matter  treated 
in  Section  V. 

In  order  to  complete  the  calculation  of  the  response  function,  even 
with  the  model  of  uniform  combustion  and  quasi-static  behavior  of  the 
gas  phase,  some  further  approximations  are  required.  There  are  three 
principal  ones  used  in  Ref.  3  and  discussed  here  in  Section  IV-A: 

1.  Combustion  in  the  gas  phase  begins  immediately  at  the  solid- 
gas  interface 

2.  The  flame  zone  is  relatively  thick,  in  a  sense  made  precise 
in  Sect  u  IV-A 

3.  Hie  fluctuations  of  energy  release  in  the  gas  phase  are  uni¬ 
form  and  hence  always  in  phase  with  the  pressure 

It  is  generally  true  that  the  flame  thickness,  Xf,  measured  from  the 
interface  to  the  outer  boundary,  can  be  assumed  to  be  large,  in  the  sense 
that  exp(mCpXf/kg)  is  very  much  greater  than  unity.  This  implies  (Eq.  58 
and  also  remarks* following  Eq.  164)  that  a  relatively  small  fraction  of 
the  heat  released  in  the  gas  phase  is  transferred  to  the  solid  by  heat 
conduction.  It  also  implies,  as  remarks  following  Eq.  58  clarify,  that 
if  in  addition  one  assumes  that  combustion  begins  at  the  surface,  the 
results  are  In  practice  restricted  to  exothermic  surface  reactions  (i.e., 
the  parameter  His  positive) .  It  appears  that  endothermic  reactions  can 
be  included  in  this  model  only  if  simultaneously  one  allows  the  combus¬ 
tion  zone  to  be  displaced  from  the  surface. 

In  Section  IV-C,  a  displaced  combustion  zone  is  treated,  although 
an  approximation  used  in  Eq.  78  implies  that  the  results  are  still  valid 
only  for  H  >  0.  Nevertheless,  the  results  indicate  that  significant 
numerical  differences  exist  between  the  calculations  of  Ref.  2  and  those 
for  the  more  realistic  case  of  a  displaced  combustion  zone.  The  dis¬ 
tinction  is  interpreted  by  means  of  the  parameter  B  in  Eq.  85.  Use  of 
the  results  of  Ref.  2  in  the  interpretation  of  data  can  therefore  be 
done  only  with  (possibly  serious)  reservations. 

The  fourth  approximation  listed  above  is  examined  briefly  in 
Section  IV-E.  It  is  supposed  that  part  of  the  fluctuation  of  reaction 
rate  is  in  phase  with  the  flame  temperature.  The  final  response  func¬ 
tion  appears  as  Eq.  106,  and  once  again  the  results  can  be  interpreted 
in  terms  of  H,  Eq.  111.  This  crude  approach  indicates  that  the  influ¬ 
ence  of  temperature  dependence  of  the  reaction  '’ate  is  less  significant 
then  that  of  a  displaced  combustion  zone. 


3 


'  '*»25»vs«wv  »x  - 


NWC  TP  4668 


All  of  these  result*  for  the  response  function  necessarily  have 
the  same  fore,  Eq,  71.  fiut  analysis  of  data  in  Ref.  7  indicates  that 
this  fona— irrespective  of  the  particular  formula  used  for  B — does  not 
accommodate  all  data.  Bence,  it  is  necessary  to  alter  one  or  more  of 
the  basic  assumptions  enumerated  earlier.  In  Section  V,  the  assumption 
of  quasi-static  behavior  is  relaxed — weakly.  An  expansion  in  the  rele¬ 
vant  frequency  parameter  (He,  Eq.  122)  is  carried  out  to  first  order,  so 
the  results  are  restricted  to  relatively  low  frequencies.  For  simpli¬ 
city,  the  calculations  are  based  on  the  model  of  uniform  combustion 
extending  from  the  solid  surface.  With  a  number  of  fairly  reasonable 
numerical  approximations,  the  response  function  eventually  takes  the 
form  shown  in  Eq.  165.  The  major  correction  associated  with  nonquasi¬ 
static  behavior  appears  in  the  imaginary  part  of  the  response  function 
and  can  be  numerically  significant.  This  conclusion  appears  to  be  con¬ 
sistent  with  observations,  but  satisfactory  proof  will  be  found  only  by 
uae  of  Eq.  165  to  analyse  data  in  a  manner  such  as  that  used,  for 
example,  in  the  work  of  Ref.  7. 

There  evidently  remain  a  number  of  significant  problems  in  respect 
to  the  unsteady  behavior  of  a  burning  solid*  Modeling  of  the  combustion 
rone  in  the  gas  phase  is  still  in  a  very  primitive  state,  and  this  is 
probably  the  severest  restriction  cn  all  the  special  cases  treated. 
Moreover,  even  for  the  simple  model  used  as  a  basis  here,  there  are 
several  permutations  not  yet  considered:  for  example,  nonlinear  behavior 
is  en  important  aspect  which  has  not  been  covered  here;  the  effect  of  a 
displaced  combustion  zone  on  nonquasi-static  and  of  course  on  nonlinear 
results  has  net  been  treated;  and  the  combination  of  nor.quasi-static  end 
nonlinear  behavior  is  bound  to  be  of  interest,  particularly  for  propel¬ 
lants  exhibiting  thick  combustion  zones.  In  any  case,  the  possibilities 
for  obtaining  analytical  results  relatively  easily  should  net  be  too 
quickly  passed  over  in  favor  of  specific  numerical  calculations. 


II.  SOLID  PHASE  AND  INTERFACIAL  MATCHING  CONDITIONS 


Since  only  the  gas  phase  td.ll  differ  from  most  earlier  works,  the 
other  pieces  of  the  problem  will  be  covered  first  and  but  briefly;  ex¬ 
panded  presentations  may  be  found  in  Ref.  1  and  8  and  sources  cited 
there.  The  model  used  is  sketched  in  Fig.  1,  with  coordinate  system 
chosen  so  that  the  origin  is  fixed  always  to  the  mean  position  of  the 
burning  surface;  the  latter  fluctuates  with  amplitude  xs.  In  this  pic¬ 
ture,  then,  cold  solid  flows  in  from  the  left  with  speed  equal  to  7, 2 
the  average  linear  burning  rate,  and  the  unsteady  heat  conduction 
(energy)  equation  for  the  solid  is 


3T 

3t 


_  3T 
pc  rc  ta 


(1) 


For  the  case  of  small  amplitude  motions ,_the  temperature  is  split  into 
its  average  value  (Independent  of  time)  T,  and  the  fluctuation  T*; 

T  »  T(x)  +  T'(x,t).  The  corresponding  solutions  of  Eq.  1  are 


(T-T  )/(f  -T  )  -  e5c 
c  sc 


(2) 


T' 


T* 

o 


(3) 


where  Tc  is  the  temperature  of  the  cold  propellant,  Ts  is  the  average 
surface  temperature,  is  the  temperature  fluctuation  evaluated  at 
x  «  £c  *  0,  and  £c  -  ‘Bc/tCg,  icc  being  the  thermal  dlffusivity,  icc  * 
kc/pcc.  The  complex  function,  X,  of  dimensionless  frequency  ft  « 
Kgtt/r2  has  real  and  imaginary  parts  given  by 


Ar  -  h{l  +  [(1  +  16 ft2)*5  +  l]*5} 


(4) 


2  This  definition  of  the  coordinate  system  will  be  altered  slightly 
in  Section  V-B. 
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xi "  iH Ul  +  16a2)i*  ’ l]k  (5) 

It  trill  be  necessary  later  to  match  the  solutions  for  the  solid  and 
gas  rhases  at  the  interface,  and  in  particular  the  fluctuations  of  tem¬ 
perature  and  heat  transfer  at  x  *  x8  must  be  computed.  For  small  potions 
of  the  surface,  these  can  be  expanded  in  Taylor  series  about  the  origin, 
and  to  first  order 


where  (  )0±,  (  )g*  refer  to  the  gas  (+)  or  solid  (-)  side  of  the  origin 
(o)  or  interface  (s).  By  use  of  Eq.  2,  3,  6,  and  7,  one  can  then  find 
the  formula  for  the  fluctuation  of  heat  transfer  from  the  interface  to 
the  solid, 


after  Eq.  13,  deduced  below,  is  inserted  for  x8.  As  shown  later  in 
Section  VI  the  same  result  can  be  found  in  a  different  way. 

There  ere  two  important  matching  conditions  which  can  be  deduced  by 
applying  conservation  of  mass  and  energy  to  a  small  control  volume  in¬ 
cluding  the  interface: 


p  (r  -  x  )  -  p  (u  -  x  ) 

C  3  gS  gS  S 


(9) 


CF  -  x  )h  _  +  (1c  4^ 


8  8- 


(s  i?J_. '  V".  -  VV + 


(>.  ®L 


(10) 


The  steady-state  form  of  Eq.  9  is  simply  an  identity,  while  that  for 
Eq.  10,  after  use  of  Eq.  2,  is 

(k,  it ' w  !c  <*.  *  V  +  L.! 


(11) 
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vh-re  Ls  -  h8+  -  he.  is  the  latent  heat  for  the  surface  reaction,  L8 
being  positive  for  an  endothermic  reaction. 

The  fluctuating  part  of  Eq.  9  gives  a  formula  for  the  surface  oscil¬ 
lations: 


a  1  -  p  /p  W 

gs  Kc 

or,  for  harmonic  oscillations,  x8  -  eiwt:,  with  to  +  0,  and  p^g  pc  «  1, 
p  x  ,  m  * 

.  _  L.  (12 

m  iu  I 

The  unsteady  part  of  the  energy  balance  in  2q.  10  leads  to 


Ls  +  5  (Cp-C)  ts* 


and  with  Eq.  8, 


/  c 

*  T  -T  L 

m  ' 

■  me 

XT  '  +  (-£  - 

4  T  ’  {  c  +  ®} 

3 

s  Vc 

> 

1  s  1  X  c  1 

m 

Finally,  an  assumption  must  be  made  in  respect  to  the  rate  of  con¬ 
version  of  solid  to  gas  at  the  interface.  It  is  common  practice  to  use 
a  pyrolysis  law  of  the  Arrhenius  form 

n_  a  -E  /RT 

n  3  v  S  8  0  8  r,,v 

“s  “  V  TS  *  (16) 


for  which  the  fluctuation  in  mass  flux  ls  given  by 

m  *  T  *  , 

■  E  =r—  +  m  2— 

Iff  T  8  p 

8 


where  E  »  a  +  E  /R  T  .  Combination  of  Eq.  15  and  17  yields 

S  8  0S 

,  /  .\  T  '  /c 


".ty'-Wt 


7 


where  A  «  (1-TC/T-)  (aa  +  E 
introduction  and  L  ■  T8/cTs 
tha  solid  and  interface  reg 


! 

* 
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III.  STEADY-STATE  SOLUTION  FOR  THE  GAS  PHASE 
WITH  UNIFORMLY  DISTRIBUTED  COMBUSTION 


The  solution  to  the  problem  was  summarized  in  Ref.  3,  but  not  in  a 
form  immediately  suitable  to  the  analysis  of  unsteady  behavior.  Thus, 
a  more  complete  discussion  is  necessary  here.  As  formulated  in  both 
Ref.  2  and  3,  the  analysis  of  the  burning  solid  is  carried  out  within 
a  wholly  "thermal  theory" — the  influence  of  diffusion  appears  nowhere 
explicitly.  The  conditions  under  which  this  Idealized  licit  is  attained 
are  discussed  in  Ref.  8  and  in  the  literature  of  flame  theory;  the  pres¬ 
sure  may  be  taken  to  be  uniform  in  space  and  the  problem  then  consists 
in  solution  of  the  energy  equation  for  the  gas  phase 


(19) 


where,  to  conform  with  the  notation  of  Ref.  2,  Qf  is  the  heat  release 
per  unit  mass  and  e  is  the  reaction  rate  (sec*’1).  For  convenience,  and 
also  in  accord  with  Ref.  2,  k  and  cp  will  be  taken  to  be  constant; 
however  (contrary  to  Ref.  2),  cp  jt  c.  The  calculations  to  be  carried 
out  here  are  based  on  the  case  when  pge  ■  w  is  constant  for  x^  x  xj 
and  zero  elsewhere,  as  shown  in  Fig.  I. 

At  the  downstream  side  of  the  flame  (x  >^Xf),  the  boundary  con¬ 
ditions  are 


T  ■  T 


f 


(x  >_  xf) 


£  -  0  (20a, b) 

and  at  the  surface,  Eq.  11  holds:  note  that  in  this  section  T  and  T  are 
Identical,  as  are  m  and  m  since  only  strict  steady  conditions  are  con¬ 
sidered.  It  should  be  noted  that  Ls  here  is  -Qs  of  Ref.  2,  and  x  here 
is  -x  of  Ref.  2.  For  the  purposes  of  this  section,  then,  the  boundary 
condition  Eq.  11  can  be  written 

(k,  ■ "  [c  «.  -  V  ’  (21) 


9 


«  **■?-.-* 
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FIG.  1.  Sketch  of  the  Model  Used  for  a  Solid  Pro¬ 
pellant  Burning  With  Uniform  Combustion  in  the  Gas 
Phase. 
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Energy  conservation  applied  to  the  gas  phase  only  gives 

(k«  f L  ■  “  'V  -  V 

Thus,  Eq.  21  and  22  lead  to 


(22) 


m  (h  .  -  hf)  ■  m  [c  (T  -  T  )  -  Q  ] 
s+  f  s  c  xs 


(23) 


Now,  by  definition  of  the  enthalpy  function. 


h 


f 


8+ 


fT 

s 

c 

ps 

h 

r 


dT 


where  hfr,  hgr  are  the  enthalpies  of  formation  at  the  reference  condi¬ 
tions  Tr  and  some  pressure  which  may  be  taken  to  be  the  pressure  at 
which  the  burning  occurs.  For  constant  specific  heats,  and  c  *  ■  c_g, 
the  definition  therefore  yields  P 


hf  -  h*e  *  -  if  +  %  <Tf  -  V  <24> 

in  which  the  heat  release  Qf  in  the  gas  phase  is  by  definition  Hsr  -  h£r 
and  is  the  as  Qf  in  t-h*  energy  equation,  Eq.  19.  The  overall 
energy  balance  for  the  gas  phase  can  be  written,  by  combining  Eq.  22  and 
24,  as 


(kg£)st'"[<if-cp<Tf-v)  «5> 

The  last  equation  is,  of  course,  simply  an  expression  for  the  over¬ 
all  energy  balance  of  the  process;  a  simple  limiting  form  arises  when 
Cp  *  c,  so  that  when  equated,  Eq.  21  and  25  give 


c 

P 


V 


Qf + 


(26) 


which  is  practically  obvious,  since  Qf  +  Qg  is  the  net  heat  release  due 
to  chemical  reactions.  But  the  point  of  this  exercise  is  to  compare 
Eq.  25  with  the  corresponding  expression  deduced  directly  from  the 
energy  Eq.  19.  Integrate  that  equation  from  the  surface  to  x  +  « 
downstream  of  the  flame  and  apply  the  boundary  conditions  Eq.  20a, b: 


11 
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(27) 


Equation  27  la  consistent  with  Eq.  25  only  if  the  reaction  rate  satisfies 
the  normalisation  condition 


e  dx  *  o 


'  o 


(28) 


In  Ref.  3  it  was  shown  that  this  condition  is  redundant  for  the  case  of 
a  very  thin  flame  front  (combustion  concentrated  at  a  plane) — it  gives 
no  information  that  cannot  be  obtained  by  other  means.  This  is  also 
true  for  the  case  of  uniform  combustion.  However,  it  is  a  convenient 
relation  to  have  at  hand,  particularly  for  the  problem  of  quasi-static 
burning  in  response  to  pressure  changes.  It  is  worthwhile  noting  that 
Eq.  28  also  leads  directly  to  a  relation  which  was  introduced  on  an 
approximation  in  Ref.  2  but  which  is  in  fact  exact.  For  suppose  that 
the  combustion  occurs  only  in  the  region  <_  x  £  Xf  and  Is  uniform 
there;  then  Eq.  28  gives 

Pg  e  (xf  -  xt)  «  m  (29) 

In  the  extreme  case  x^  ■  0, 

W  =  P_  E  -  (30) 

8  *f 

a  relationship  which  was  adopted  as  an  approximation  just  before  Eq.  8 
of  Ref.  2. 


As  a  convenience3  for  integrating  Eq.  19,  it  is  helpful  to  define 
a  new  variable  C  (called  £  in  Ref.  3): 


C 


me 

k. 
e  8 


dx 


so  that 


(31) 


d_  d_  d 

dx  dx  d£  k  5  dC 

8 


(32) 


3  And  only  a  convenience — the  steps  leading  to  Eq.  32  are  obviously 
unnecessary,  but  using  £  Instead  of  x  does  simplify  writing  subsequent 
manipulations . 
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With  this  transformation,  the  left-hand  side  of  Eq.  19  is 


me 


—iril _ 2  t  —  (me  C—  = _ 2 —  c  ±i 

kg  5  dC  ks  C  dC  lOCp  C  dCi  kg  ?  d?2 


g 


and  the  whole  equation  can  be  written 

-=23F-*2 

with  t  *  T/'fs,  and  the  eigenvalue  A2  is  defined  as 


(33) 


A2 


mzc  ZT 


(34) 


Since  w  ■  Pgl  is  here  taken  to  be  nonzero  only  in  the  interval 
£  x  <  r,^  where  it  is  constant.  A2  is  also  nonzero  and  constant  only 
in  the  same  region.  Thus,  Eq.  33  is  easily  integrated  to  give  dr/dc  and 
t  as  follows: 


dr 

dC 


C1 


C1  C 


+  c. 


Cs  i  i  Ci 


(35) 

(36) 


dr 

d  C 


c 


3 


t  *  A2  tn  C  +  c^  ?  +  c^ 


(37) 

(38) 


Cf  1  Cf  (39a, b) 


There  are  eight  conditions  to  be  satisfied:  the  boundary  values 
T  «  Ts  and  Eq.  21  at  the  surface,  Eq.  20a, b  downstream  of  the  flame 
(x  ■+■  *) ,  and  continuity  of  temperature  and  heat  transfer  at  the  bound¬ 
aries  of  the  flame,  x  *  and  x  *  xj.  As  part  of  the  given  information 
defining  the  problem,  w  and  an  "ignition  temperature,"  Tif  at  which 


dr 

de  “  C5 


r  «  c5  C  +  c6 
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combustion  begins,  muse  be  specified.  From  Tj,  x^,  or  Ci  can  be  calcu¬ 
lated  from  Eq.  36.  There  ere  therefore  eight  quantities  determined 
by  the  eight  conditions  listed  above:  the  six  constants  of  integration, 
Cj,  . ...  eg,  the  flame  thickness  Xf,  and  finally  the  eigenvalue  A2.  The 
end  result  is  an  expression  for  A2  which  can  be  solved  to  give  the  mass 
flux  or  linnar  burning  rate  as  a  function  of  pressure  and  other  vari¬ 
ables  in  the  problem. 

A  direct  way  to  apply  the  various  conditions  is  to  find  the  temper¬ 
ature  profiles  for  the  upstream  and  downstream  slde3  of  x  -  x*  and  match 
there.  First,  the  surface  condition  on  heal  transfer,  Eq.  21,  is  in  the 
normalized  variables: 


«.  (fL  •  rr  [c(Ts  -  V  -  o. 

\  Vs*  p  s 


The  correct  values  are  taken  by  the  profile  Eq.  35-36  at  the  surface  if 


fdt\ 

'  W, 


c2  *  %  - 


/  dr\ 

w. 


so  that  between  the  Interface  and  the  upstream  edge  of  the  flame, 


t  ■  t  + 
8 


« -  v  (If). 


(c8  <  ?  <.  ?) 


where  (d-r/d?)8+  is  given  by  Eq.  40. 

The  conditions  Eq.  20a, b  far  downstream  require 


Cj  «  0 


£  -  ■  Y 

c6  Tf 


(43a, b) 


Hence,  continuity  of  temperature  and  heat  transfer  at  ( 
satisfied  if  in  Eq.  38  and  39 


'3 


c^  -  rf  +  A2  (1  -  In  Cf) 


(44a, b) 
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and  the  temperature  profile  in  the  flame  zone  is 


t  ■  T, 


-4(r)-Rj 


<c±  <  c  <.  cf> 


(45) 


which  must  be  matched  to  Eq.  42  at  (  »  Continuity  of  temperature 

there  gives 


Tf  "  Ti _ 

£n  it  _  (i  _  Si) 

Ci  cf 


(46) 


and  continuity  of  heat  transfer  Implies 


Cf  Ci  (dt\ 

rf  *  ?i  W 


(47) 


8+ 


It  is  easily  verified,  when  due  account  is  taken  of  the  definitions 
involved,  that  Eq.  47  is  exactly  Eq.  7  of  Ref.  2. 

For  completeness,  it  may  be  noted  that  the  normalization  condition, 
Eq.  29,  in  dimensionless  variables  is 


(Qf/c  T  ) 


(48) 


Any  two  of  Eq.  46-48  give  an  equation  which  must  be  solved  to  find 


SMIL-© 


(49) 


(Observe  for  use  in  Eq.  46  that  according  to  Eq.  25,  Tf  “  *  'C~Y" 


-  (dt/dO^O 


p  s 


Numerical  calculation  of  the  steady-state  problem  proceeds  as 
fellows: 


(1)  The  quantities  Tf,  Qs,  Qf  are  assumed  known  by  thermodynamics 
or  otherwise,  while  the  temperature  of  the  cold  propellant  Tc  can  be 
chosen  as  desired.  The  reaction  rate  w  must  be  specified  as  a  function, 
say,  of  pressure  and  average  temperature. 

(2)  Choose  a  value  for  Ts;  consider  the  case  ng  ■  0  so  that  m  can 
be  computed  from  the  pyrolysis  law,  Eq.  16 


15 


u  : 

.v  <  • 


W~W  'i? 


ttLgiftM 


_  «.  "W. 

*"B.  T.  e 


wh«r«  Bg,  4|,  E#  are  assumed  known. 

(3)  In  the  purely  steady-state  problem,  the  Interface  can  be  fixed 
at  the  origin,  so  x,  *  0  and  5a  »  1;  (dt/d?)^  can  then  be  computed  from 
Eq.  40  with  c_  -  1: 


(s?L  -  rr  >C<W  - «. 

\  /8+  p  S 


(50) 


(4)  The  position  of  the  inner  edge  of  the  flame  follows  from  the 
specified  value  of  used  in  Eq.  42,  again  with  ?s  •  1: 


r  m  ,  +  VIl 

Zi  (dr/dc)  . 


(51) 


(5)  The  position  of  the  outer  edge  of  the  flame  can  be  found  from 
the  transcendental  Eq.  49. 

(6)  The  value  of  the  eigenvalue  can  then  be  calculated  from,  say, 
Eq.  48: 

P  * 

from  which,  since  m  was  already  calculated  in  step  (2),  a  number  is 
found  for  w.  Since  the  reaction  rate  is  generally  a  function  of  pres¬ 
sure  (specified),  this  last  step  gives  the  pressure. 


Thus,  the  net  result  is  the  variation  of  mass  flux  or  linear  burn¬ 
ing  rate  with  pressure.  If  n8  i  0  in  Eq.  16,  then  an  iterative  calcu¬ 
lation  is  required:  a  value  of  pressure  must  be  assumed  initially  to 
compute  m  from  the  pyrolysis  law.  The  value  computed  in  step  (6)  must 
then  be  compared  with  the  value  assumed  in  step  (1)  and  the  calculations 
repeated  until  satisfactory  agreement  is  achieved. 

This  is  an  interesting  problem  to  complete,  for  comparison  with 
experimental  data  for  the  burning  of  composite  propellants:  a  thorough 
discussion,  with  numerical  results,  is  being  prepared.  The  point  of 
Including  the  discussion  here  is  to  clarify  just  what  must  be  known  to 
solve  the  problem,  and  just  what  can  be  calculated,  for  a  clear  under¬ 
standing  of  steady-state  burning  in  this  respect  is  helpful  in  analysis 
of  the  quasi-static  behavior. 
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IV.  QUASI-STATIC  BEHAVIOR  OF  THE  GAS  PHASE 


When  the  pressure  varies  In  time,  but  is  still  assumed  to  be  uni¬ 
form  in  space,  the  equations  of  motion  for  the  gas  phase  are 


3p  .  3m 
9t  3x 


0 


(52) 


3T  .  3T 
PCp  3g  mCp  3x 


+  QfPg* 


(53) 


It  is  a  formal  implication  of  the  physical  assumption  of  quasi-static 
behavior  that  the  time  derivatives  in  the  differential  equations  of 
motion  can  be  neglected.  This  is  based  on  the  idea  (perhaps  supported, 
but  not  necessarily  wholly  justified,  by  order  of  magnitude  estimates) 
that  the  processes  within  the  gas  phase  respond  essentially  instan¬ 
taneously  to  changes  of  pressure  and  of  boundary  values — temperature, 
heat  transfer,  and  positions  of  the  boundaries. 

The  origin  of  the  time  dependence  is  of  course  the  pressure  field 
external  to  the  flame,  but  the  assumption  of  instantaneous  response 
does  not  mean  that  other  variables  in  the  problem  change  in  phase  with 
the  pressure.  In  fact,  a  significant  lag  is  introduced  by  the  behavior 
of  the  thermal  wave  in  the  solid;  this  Influences  the  behavior  of  the 
gas  phase  mainly  through  the  heat  transfer  from  the  interface  to  the 
solid,  expressed  in  Eq.  8.  But  in  consequence,  the  fluctuations  of 
surface  temperature,  flame  temperature,  flame  thickness,  and  reaction 
rate  also  cannot  generally  take  place  in  phase  with  the  pressure  except 
in  the  limit  u  ■  0.  Nevertheless,  the  solution  to  the  strictly  steady- 
state  solution  can  still  be  used  if  the  influence  of  fluctuations  is 
properly  accounted  for. 

An  example  may  serve  to  clarify  what  is  meant  by  "quasi-static" 
when  time  lags  are  present.  Consider  the  simple  case  of  heat  conduc¬ 
tion  between  two  planes  at  different  temperatures  T^,  T2  and  separated 
by  a  distance  L.  The  temperature  profile  in  the  steady  state  is  linear, 

T(x)  -  TX  +  (Tj-Tj)  ~  (5A) 

where  the  temperature  T  ■  Tj  at  x  «  0.  The  assumption  of  quasi-static 
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behavior  when  small  fluctuation  Tj.*,  Tj'*  L’  occur  Implies  that  the 
fluctuation  of  the  local  temperature  profile  can  be  calculated  simply 
by  computing  Its  perturbation  according  to  the  steady-state  solution 
Eq.  54.  Write  T  •  T  +  T',  Tj  -  +  Tj/,  etc.,  and  since  the  barred 

quantities  obey  Eq.  54,  one  finds  easily  a  formula  for  T': 

T’(x)  -  Tj*  +  (T2'  -  T1*)  i  -  (?2  -  Tjl)  |  (=H  (55) 

The  total  temperature  profile  T  «  T  +  T*  as  well  as  T  and  T*  separately 
are  still  linear  functions  of  the  distance  from  the  colder  (T^)  plate 
but  the  ends  of  the  profile  need  not  jiggle  back  and^forth  in  phase. 

It  is  entirely  permissible,  for  example,  that  Tj_'  ■  T^'  cos  wt  and  ?2'  * 
fo'  cos  (ut-4)  so  that  the  fluctuation  in  temperature  of  thejiotter 
plate  lags  in  time  by  the  constant  angle  4*  Moreover,  fj/,  T2 '  could  be 
specified  as  functions  of  the  frequency,  to.  Similar  statements  apply 
to  the  variable  separation  of  the  planes.  The  point  is  that  the  steady- 
state  solution  fixes  the  functional  form  of  the  solution — in  this  case 
linear — while  phase  lags  and  dependence  on  frequency  can  be  introduced 
through  the  boundary  values.  A  truly  nonsteady  calculation  would  of 
course  lead  to  a  thermal  wave  in  the  region  between  the  planes.  The 
distinction  between  these  cases  is  illustrated  in  Fig.  2. 
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FIG.  2.  An  Illustration  of  the  Difference  Between  Quasi-Static 
and  Truly  Unsteady  Temperature  Profiles  in  a  Solid  Slab. 
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The  steady-state  solution  for  the  gas  phase  found  in  the  preceding 
section  can  be  handled  in  exactly  the  same  way.  Although  the  form  of 
the  profile  is  fixed  by  Eq.  42  and  45  within  the  gas  phase,  the  effect 
of  the  thermal  wave  in  the  solid  is  to  force  both  a  phase  lag  and  a  de¬ 
pendence  on  frequency  for  practically  all  variables  in  the  problem.  In 
particular,  and  this  point  will  be  mentioned  again,  the  fluctuation  of 
the  flame  thickness  Xf  forces  a  fluctuation  in  the  ratio  w/m  *  1/xf 
vt.ich  is  not  in  phase  with  the  pressure. 

It  is  perhaps  best  to  begin  the  following  manipulations  by  summar- 
iz  the  results  already  obtained  and  required  here.  What  is  needed  is 
a  :  i aula  for  the  fluctuation  of  heat  transfer  from  the  gas  phase  to  the 
int  rface,  for  use  in  the  energy  matching  condition,  Eq.  18.  If 
(kg  oT/dx)s+  is  expressed  as  a  linear  combination  of  contributions  due 
to  ms',  Ts',  and  p',  then  with  Eq.  17,  there  is  sufficient  information 
to  give  a  formula  for  the  response  ratio  mg'/p'*  Thus,  one  already  has 


(17) 


(18) 


For  easier  comparison  with  Ref.  2,  L  *  L  /cT  has  been  replaced  by 

s  s 

H  -  Q  /c  (T  -  T  )  «  -  L/(l  -  T/T) . 
s  s  c  c  s 


When  the  assumption  of  quasi-static  behavior  is  made,  there  are  two 
ways  of  arriving  at  a  formula  for  kg(3T' /3x)g+,  both  giving  the  same 
result.  Onejway  is  to  linearize  the  two  equations,  Eq.  52  and  53,  by 
writing  T  *  T  +  T'  and  solving  for  the  fluctuations  with  the  time  deriva¬ 
tives  ignored.  This,  as  shown  later  in  Section  V,  is  really  the  first 
term  of  an  expansion  in  powers  of  a  frequency  parameter.  The  second  way 
is  simply  to  linearize  the  steady-state  solution,  as  done  above  for  heat 
flow  between  two  plates.  The  equivalence  of  these  two  methods  is  readily 
demonstrated  and  is  shown  below. 


A.  Quasi-Static  Results  from  the  Steady-State  Solution 

This  means  of  finding  the  quasi-static  formulas  has  been  widely 
used,  sometimes  incorrectly:  any  steady-state  solution  for  the  gas  phase 
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can  be  used  as  a  basis.  The  steady-state  problem  of  Section  III  from 
which  the  quasi-static  formula  for  (k_  3T'/3x)g+  may  be  found  has  yielded 
tha  following :  8 

(1)  The  temperature  profile 


■ +  (t  -  '•>  (iuL  - 5  -  'i1 
\  /  0+ 


T  -  T, 


-  A2  *n  f  -  1  -  f-  (?i  <  ?  i  Cf> 


2  2=: 


(2)  Three  expressions  for  the  eigenvalue  A2  =  ka  wQf/m  cD  t 


Tf  “  T1 


Vl  (dr) 
'f  -  C1 


(3)  A  transcendental  equation  for  the  flame  thickness,  found  by 
equating  Eq.  47  and  48 


?■>*(  *L--(?) 


(4)  The  overall  energy  balance  for  the  gas  phase  Eq.  25 


k8  (£)  ^  -  -  -  <=„  <Tf  -v1 


which  in  dimensionless  variables  is 


(dr  \  . 

\«L 


<Tf  -  T.) 


(25a) 
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(5)  A  normalization  conditional  Eq.  29,  on  the  reaction  rate  w, 
which  la  implied  by  Eq.  25  and  which,  rewritten,  la  Eq.  48: 


(29) 


Not  all  of  the  above  are  independent  re8ults:  by  u8e  of  the  expressions 
for  the  temperature  profile  and  the  energy  balance  Eq.  25a,  any  two  of 
the  formulas  for  the  eigenvalue  can  be  derived  from  the  third.  The  model 
has  been  solved  exactly,  but  numerical  results  cannot  be  obtained  until 
the  reaction  rate  w  is  independently  specified. 

The  fluctuation  of  heat  transfer  is  found  by  linearizing  in  the 
usual  way:  T  «  T  +  T',  m  ■  m  +  m',  etc.  A  question  arises  whether  the 
values  indicated  as  (  )g+  here  are  really  those  at  the  Interface  or 
whether,  as  for  the  solution,  Eq.  3,  to  the  solid  phase,  they  represent 
values  on  x  •  0,  from  which  the  correct  values  at  the  interface  must  be 
found  by  using  the  linearized  boundary  conditions  Eq.  6  and  7.  The 
answer  is  that  within  the  quasi-static  assumption,  the  terms  involving 
xs  must  be  ignored.  The  reason  for  this  is  that  in  the  limit  <n  *  0, 
the  fluctuating  part,  x8,  of  the  surface  motion  is  zero,  so  that  x8  is 
a  constant  which,  by  free  choice  of  the  origin,  may  be  taken  to  be  zero. 
In  other  words,  it  is  inappropriate  to  use  the  formula  x8  «  iux8  and 
hence  in  the  limit  w  ■  0,  x8  cannot  be  found  from  Eq.  13.  Incidentally, 
it  is  easily  verified  that  the  correct  quasi-static  limit  (i.e.,  o»  *  0) 
does  follow  from,  for  example,  Eq.  18  for  the  solid  phase,  even  though 
that  relation  involves  use  of  Eq.  13.  The  reason  for  this  is  that  the 
material  of.  the  solid  phase  is  assumed  to  move  always  at  the  uniform 
speed  r  ■  m/pc.  In  the  treatment  of  the  gas  phase,  fluctuations  of  the 
speed  of  the  gas  are  necessarily  included,  and  hence  neglect  of  the 
terms  in  x8  must  be  enforced  as  a  separate  statement,  part  of  the  quasi¬ 
static  analysis.  A  more  careful  handling  of  the  surface  boundary  con¬ 
ditions  is  requiipd  if  the  assumption  of  quasi-static  behavior  is  re¬ 
laxed,  or  if  nonharmonic  changes  of  pressure  are  considered.  The  point 
here  is  that  so  far  as  periodic  quasi-static  variations  are  concerned, 
values  denoted  by  (  )„+  are  Identical  with  values  (  More  detailed 

considerations  appear  in  Section  V. 


To  find  an  expression  for  (dT'/dx)^,  consider  Eq.  47,  with  the 
definition  of  A2: 


Vkl 

2  2_ 
m  c  T 
P  8 


which  gives,  since  ( 

8 


1  for  x  ■  0, 
s 
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A 8  already  renarked,  this  can  be  found  by  direct  integration  of  the 
energy  equation  and  appears  as  Eq.  7  of  Ref.  2  and  Eq.  87  below.  Before 
considering  refinements,  the  simple  esse  of  uniform  combustion  beginning 
at  the  surface  will  be  discussed,  the  problem  handled  in  Ref.  2;  thus, 

Z±  a  1.  Moreover,  the  approximation  that  Cf  >>  1  will  be  introduced  at 
this  point,  and  used  in  almost  all  the  remaining  discussion.  Hence, 

Eq.  56  becoaes 


(57) 


A  useful  interpretation  of  "large  £f"  follows  from  Eq.  47  and  48, 
which  show  that  for  large  ?£, 


and  hence  the  heat  transfer  to  the  solid  from  the  gas  is  small  compared 
to  the  energy  transferred  by  convection.  Moreover,  with  Eq.  40,  this 
relation  becomes 


(58) 


where  H  *  Qa/c  (T8-Tc)  will  appear  as  an  Important  parameterjLn  the  final 
results  for  the  response  function.  For  practical  cases,  Tc/Ts  -  1/3,  so 
that  one  must  choose  H  near  unity  to  be  consistent  with  the  approximation 
that  Hr-  is  large.  However,  this  is  not  especially  troublesome,  providing 
H  >  0,  since  Cf  can  be  considered  large  when  £n4  £  2.  Then  since 
Cp/c  -  lA  1-R  -  Cl/2) / (2/3)  *  0.75  and  H  »  0.25;  in  other  words,  the  approx¬ 
imation  Cf  »  1  accommodates  positive  values  of  H  in  the  range  0.25  -  1. 


A  rough  numerical  estimate  indicates  that  may  indeed  be  quite 
large.  Since  »Cp  x£/kg  «  (pccp/kg)  rxf  *  (cpkc/ckg)  rxf/xc,  and 

cpkc/ckg  *  1,  one  has  -  exp  (rxj/xc) .  A  representative  value  of  kc 

is  10”  ^  cm^/sec  and  for  r  ■  1  cm/sec,  -  exp  (1000  Xf)  with  xj  in 

centimeters.  ^JSven  for  Xf  as  unrealistically  small  as  20  microns  *  0.002 
centimeters,  Xf  * 


Vi 
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If,  on  the  other  hand,  H  is  negative,  then  instead  of  Eq.  58  one 


i-a 


(1  +  (-H))  - 


+r4i~hr) 

tet/z±) '  «±  cf  / 


(58a) 


and  it  is  apparent  that  H  <  0  and  large  values  of  £f  are  not  compatible 
unless  (£  is  also  fairly  large.  The  treatment  of  Ref.  2  is  therefore 
restricted  to  positive  values  of  H  (i.e.,  exothermic  surface  reactions) 
because  the  assumptions  are  made  that  £f  is  large  and  that  **  1  (i.e., 
combustion  begins  at  the  surface) . 

The  quasi-static  formula  for  the  fluctuation  of  heat  transfer  is 
found  simply  by  linearizing,  Eq.  57: 


(t  a)  2ti!( 

V  8  y  ;  / 

-;«i,  a2  (^-  -  al) 

p  8  \w  m  l 


where,  as  usual.  A*  «  Qf  kg  w/(cpm)  T  *  (Qf/c  Tg)  tfn  Cf)  *  by  Eq.  48. 
The  fluctuation  of  mass  flue  is  related  to  the  pressure  and  surface 
temperature  fluctuations  by  the  linearized  pyrolysis  law,  Eq.  17. 

A  second  formula  for  the  heat  transfer  at  the  surface  may  be  deduced 
by  linearizing  Eq.  25  for  the  overall  energy  balance  of  the  gas  phase: 


-  .-1 


ML  ■” 


'  [Q*  -  c  (T.  -  T  )]  -  me  (T  '  -  T  ') 
xf  p  f  s  p  f  s 


As  later  calculations  will  show,  this  relationship  is  used  to  compute  the 
fluctuations  of  flame  temperature,  Tf',  once  the  remainder  of  the  problem 
has  been  solved. 

The  preliminary  calculations  are  completed  by  finding  a  formula  for 
the  fluctuation,  w' ,  of  reaction  rate.  First,  following  Ref,  2,  it  may 
be  most  simply  supposed  that  w'  is  the  fluctuation  of  the  steady-state 
reaction  rate  w(p)  which  was  computed  in  the  work  in  the  following  way. 
The  average  heat  transfer  from,  the  solid  to  the  gas  is  given  by  the 
results  above  as  Eq.  58.  But  from  the  matching  condition  at  the  inter¬ 
face  and  the  temperature  distribution  in  the  solid  Eq,  11, 


Ml 


-  m  [c  (T  -  T  )  -  Q  ] 
s  c  s 
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Equating  these  expressions  gives 


a2  [c  (T  -  T  )  -  Q  ] 

sc  s 


(61) 


It  is  supposed  that  the  right-hand  side  can  be  expressed  as  a  function 
of  pressure  only.  This  is  accomplished  by  assuming  that  the  linear 
burning  rate  may  be  approximated  by  the  law  m  ■  apn  and  that  an  approxi¬ 
mation  to  the  Arrhenius  law  is 

a  -  b(Tg  -  Tc)8  p”8  (62)  *♦ 


so  that  T  -  T 
s  c 


n-n 

(p  8a/b) 


l/s. 

9 


hence. 


Eq.  61  reads 


(63) 


In  this  way  w  is  given  by  Eq.  63  as  a  function  of  pressure  only,  and  has 
been  so  constructed  as  to  satisfy  the  equations  for  steady  burning  for 
the  case  a  -  apn  and  a  pyrolysis  law  approximated  by  Eq.  62. 


The  assumption  in  now  made  that  even  during  oscillations,  the  re¬ 
action  rate  is  given  Instantaneously  by^  the  same  function  of  pressure, 
so  that  w'  is  computed  by  setting  p  **  P  +  p'  in  Eq.  63  and  retaining 
only  terms  linear  in  p'.  This  implies  that  the  local  fluctuation  of 
energy  release,  Qfw',  is  also  always  in  phase  with  the  pressure  fluctu¬ 
ations. 


Now  the  temperature  fluctuations  are  not  in  phase  with  the  pressure. 
For  example,  the  fluctuation  of  temperature  at  the  edge  of  the  flame, 

Tf ' ,  follows  from  Eq.  60: 


where  (k~  ST'/Sx)^  is  given  here  by  Eq.  59.  Thus,  if  the  local  reaction 
rate  is  at  all  sensitive  to  temperature,  there  must  be  a  part  of  w*  which 
(contrary  to  the  assumptions  used  in  Ref.  2  and  expressed  above)  is  not 
in  phase  with  pressure  oscillations.  The  result  that  w*  is  entirely  in 


**  Note  that  ng  »  0  in  Ref.  2. 
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phase  with  pressure  is  a  consequence  of  two  assumptions :  (1)  that  the 

reaction  rate  is  everywhere  in  the  gas  phase  a  function  of  pressure 
only,  and  (2)  quasi-static  behavior. 

Even  though  one  may  represent  the  steady-state  global  combustion 
process  as  a  function  of  pressure,  uniform  in  space,  this  does  not  mean 
that  locally,  within  the  gas-phase  flame,  one  can  always  simply  ignore 
other  variables.  In  the  analysis,  the  influences  of  those  variables 
are  smeared  out  in  the  steady  state  to  give  the  uniform  function  of 
pressure  w(p)  derived  above.  In  fact,  it  may  be  better  to  regard  the 
calculation  as  a  means  of  relating  the  overall  rate  of  energy  release, 
Qfwxf ,  to  the  linear  burning  rate.  Treatment  of  this  as  a  function  of 
pressure  only  amounts  to  finding  what  the  (perhaps  fictitious)  uniform 
function  of  nressure  would  have  to  be  in  order  to  match  an  observed 
burning  rate,  which  is  here  approximated  as  apn. 

The  view  that  the  reaction  rate  depends  almost  entirely  or  pressure 
is  based  on  the  idea  that  the  overall  burning  rate  is  diffusion  con¬ 
trolled.  Nevertheless,  locally  the  homogeneous  gas-phase  reactions 
must  respond  to  temperature  changes,  and  of  course  the  process  of  diffu¬ 
sion  is  itself  dependent  on  temperature.  Moreover,  under  unsteady  con¬ 
ditions,  there  are  likely  to  be  fluctuations  of  mixture  ratio.  Both  of 
these  will  affect  the  local  reactioti  rate  and  will  therefore  contribute 
pieces  to  w'  which  are  not  in  phase  with  the  pressure  change,  even 
though  the  processes  are  treated  in  a  quasi-static  manner.  Such  a 
possibility  is  not  strictly  consistent  with  the  assumption  that  the  re¬ 
action  rate  is  uniform  in  space.  That  assumption  necessarily  precludes 
consideration  of  dependence  of  w  on  variables  which  exhibit  spatial 
variations. 

Thus,  one  concludes  that  if  one  requires  the  fluctuating,  as  well 
as  the  steady-state,  reaction  rate  tc  be  uniform  in  space,  one  is  seri¬ 
ously  restricted;  the  energy  release  is  in  phase  with  the  pressure.  In 
fact,  that  conclusion  holds  even  if  the  assumption  of  quasi-static 
behavior  is  relaxed  as  shown  in  Section  V.  The  calculation  of  Ref,  2 
therefore  cannot  be  usefully  extended  without  altering  certain  approxi¬ 
mations  on  which  that  work  is  based. 

In  summary,  the  three  main  approximations  which  qualify  the  results 
of  Ref.  2  are: 

(1)  That  combustion  in  the  gas  phase  Degins  immediately  at  the 
interface 

(2)  That  the  flame  zone  is  relatively  thick,  thus,  as  the  remarks 
following  Eq.  58a  show,  excluding  consideration  of  endothermic  surface 
reactions  unless  approximation  (1)  is  relaxed 

(3)  That  the  fluctuations  of  energy  release  are  uniform  in  the  gas 
phase  and  hence  always  in  phase  with  the  pressure  fluctuations. 
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The  analysis  of  Ref.  2  is  Incomplete  so  that  the  approximation  (1)  could 
not  be  corrected;  a  procedur*  for  doing  so  is  developed  below  in  Section 
IV-C.  Correction  of  (2)  to  include  endothermic  reactions  can  be  done 
without  essential  difficulty,  but  it  is  lengthy  and  will  not  be  included 
here.  Alternatives  to  the  approximation  (3)  are  more  difficult  to  con¬ 
struct  but  may  be  more  important  than  relaxation  of  (1) .  Some  remarks 
as  this  question  are  given  in  Section  IV-D. 

If  Eq.  63  is  used,  the  fluctuation  of  reaction  rate  is 


where 


W  -  n 


c 

2(1-H)  +  -2- 
c 


(l-ng/n) 

5 

J 


(64) 


(65) 


and  H  «  Qg/c(Ts-Tc).  The  factor  cp/c  in  cp/c)(n-n8)/s  arises  from  the 
fluctuation  of  Qs,  Qs*  -  -(cp-c)Ts\  which  appeared  also  in  Eq.  14. 


B.  Calculation  of  the  Response  Function 
for  the  Simplest  Case 

Before  further  consideration  of  the  approximations  noted  above,  it 
is  useful  to  compute  the  response  function  for  the  simplest  case,  thereby 
showing  most  clearly  how  the  various  pieces  of  the  problem  fit  together. 
For  this  purpose,  all  that  is  required  is  contained  in  Eq.  17,  18,  59, 
and  64.  Substitution  of  Eq.  64  into  Eq.  59  gives  what  amounts  to  the 
net  result  of  analysis  of  the  gas  phase: 


-L_  (|l)  '.(i-L.)i,£l-f£A2aL 

mcT  '  X  /  \  T  /  p  C  m 

8  8+  S  V 


Now  substitution  of  this  formula  into  the  energy  matching  Eq.  18,  and 
use  of  the  pyrolysis  law,  Eq.  17,  to  eliminate  Ts'/Ts,  leads  to 


, (AW  +  -£  n  }  +  n  (A-I) 
iu  /m  _  \ _ c  s/  s  _ 

p*/p  A  +  y  +  EA2  -HA  +  -2.  -  1 

AC  C 


(67) 


26 


NWC  TP  4668 


Now  for  Che  assumed  form  of  the  mean  burning  rate,  m  *  ap  ,  Che  limiting 
value  of  the  above  ratio  for  w  *  0  must  be  n  (everything  changer;  in  phase 
with  the  pressure)*  Hence,  since  A  1,  this  condition  requires 


c  c 

A(l-H)  +  — ^  EA2  +  — 
c  c 


and  with  W  given  by  Eq.  65, 


AW  +  — n 

c  s 


i  ea2  +  i  (x  .  .  *  a_H>  +  i  i(l-i) 


(68) 


That  this  is  identically  satisfied  can  be  shown  in  the  following  way. 
First.-  the  index  s  (called  m(!)  in  Ref.  2)  is  in  fact  exactly  equal  to 
the  more  commonly  used  A.  To  see  this,  write  the  linearized  form  of  the 
pyrolysis  law,  Eq.  62, 


m 


1-T 


P 


By  comparison  with  Eq.  17,  obviously  s  *  E(l-Tc/Ts)  »  A.  Thus,  Eq.  68 
requires  A(l-H)  ■  cpEA2/c.  But  this  is  merely  the  mean  energy  balance 
at  the  surface,  written  in  a  slightly  obscure  form.  For  with  Eq.  57 
substituted  into  Eq.  11  with  Ls«  -Qs,  one  has 

Q,  k  w 

-  me  (T  -  T  )  (1-H) 
me  sc 

P 

and  by  use  of  the  definitions  of  A  and  A2,  multiplication  by  E  gives 


c 

EA2  -  A  O-H)  (69) 

c 


Thus,  Eq.  67  has  the  correct  limit  "built-in";  this  will  always  be  the 
case  if  correct  use  is  made  of  the  steady-state  energy  balance  (see 
Ref.  5,  Eq.  40,  for  example).  This  formula  can  be  put  into  the  standard 
form  suggested  in  Ref.  1  by  defining 


B 


n 


8 

An 


c 

2 (1-H)  + 

c 


1 

A 


(70) 


and  in  view  of  Eq.  68,  Eq.  67  can  be  written 
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.,/-  nAB  +  n  (X  -  1) 

S_il  - - 5 - 

p’/p  X  V  j  -  (1+A)  +  AB 


(71) 


It  is  clear,  therefore,  that  there  are  only  two  significant  parameters, 

A  and  B,  other  than  iCc/r2  in  X  and,  of  course,  n.  Note  also  that  it  is 
really  H  that  is  changed  to  change  B,  and  that,  by  the  steady-state  energy 
balance  at  the  surface,  Eq.  69,  1-H  can  also  be  interpreted  as  a  combin¬ 
ation  of  gaa  phase  properties: 

c  ,  Q,  k  w 

1-H  ■  J2  - - -  -*gr  (72) 

i*vt,  ‘-v  % 

The  limit  H  •  1  probably  cannot  be  reached  in  practical  cases  because 
this  implies,  according  to  Eq.  11,  that  no  heat  is  transferred  from  the 
gaa  to  the  solid — the  heating  of  thejsolijd  is  then  due  entirely  to  the 
exothermic  reaction  at  the  surface:  T8  ■  Tc  +  Qs/c.  Although  this  does 
not  violate  conservation  of  energy,  it  is  an  unlikely  dynamical  state. 

The  fact  that  Eq.  72  requires  m  -*•  •  is  not  significant  because  that  re¬ 
sult  explicitly  requires  a  nonzero  value  of  heat  transfer  from  the  gas 
to  the  solid.  Negative  values  of  H  (endothermic  surface  reactions)  are 
of  course  accommodated  by  Eq.  72. 


C.  Combustion  Zone  Not  Extending  to  the  Surface 


Since  the  assumption  that  the  gas  begins  burning  immediately  upon 
issuing  from  the  surface  appears  to  be  a  serious  restriction,  it  is  use¬ 
ful  to  examine  a  simple  way  of  avoiding  it.  Suppose,  for  simplicity, 
that  the  ratio  of  the  mean  position  of  the  inner  ed&e  to  that  of  the 
outer  edge  is  a_constant,  6  :  x^  *  Bxf .  Then  -  Cf®  and  /Z±  - 
fKCf'/Cf).  If  Cf1”®  »  1  still,  then  Eq.  57  is  replaced  by 


L  il\  .  qf  k?  1  _w_ 

\  *  Cp  *  5f* 

and  the  quasi-static  formula,  Eq.  59,  is  replaced  by 


•  me  T 
P  « 


a2  ZT-(~ 
CfB  \v 


The  argument  leading  to  Eq.  63  now  provides 


( 


w  » 


(73) 


(74) 


(75) 
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in  which  must  be  expressed  as  a  function  of  pressure  from  the  steady- 
state  problem  (see  Eq.  78  below). 


If  one  follows  the  argument  used  in  Kef.  2  and  in  Section  IV-B, 
only  the  fluctuation  w'  in  Eq.  74  should  be  computed  from  the  steady- 
state  relation  w(p)  given  by  Eq.  75.  Obserye  that  if  one  assumes  that 
the  flame  thickness  also  responds  instantaneously  and  in  phase  with  the 
pressure,  then  w/tf^  in  Eq.  73  coula  be  linearized  by  use  of  Eq.  75  and 
the  final  result  for  the  response  function  would  be  the  same  as  Eq.  71. 
Thus,  in  Eq.  74,  must  be  replaced  by  use  of  the  quasi-static  linear¬ 
ization  of  Eq.  29: 


w 


m 


xf"xi 


1 

•^n(?f/ci) 


1 

(i=KETcf 


which  gives 


so  that  Eq.  74  is 


(76) 


«  m  c  T 

p  8 


(1  +  8  £n  ?f)  ~ 
w 


-  (1  +  2  6  ta  ?f) 

m 


(77) 


Now  the  relation  for  £f(p),  the  true  steady-3tate  flame  thickness, 
is  easily  obtained  from  Eq.  49  and  50,  with  Eq.  62  used  for  Ts  -  Tc. 

One  finds,  with 


,  a  n-nnv 
(c  r  p  3) 


1/s 
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and  again  for  large  Cf  (•  circumstance  which  is  likely  always  to  be 
true) , 


Cf*  *n  ?f  * 


Qf / (1—0) 


(a  n-na\ 

hP  8) 


1/8 


(78)' 


~  Q. 


This  gives,  for  use  in  Eq.  75, 


«*' 


-  Z 


(79) 


where 


c  /c  /n-n  \  (n  (( 

iV  (t*) 


(1+8  tn  ?f) 


(80) 


The  factor  cp/c  once  again  comes  from  the  fluctuation  of  Qg.  Note  that 
Eq.  79  shesrethat  a  decrease  of  flame  thickness  accompanies  an  increase 
of  pressure,  which  is  true  for  any  value  of  8.  Linearization  of  Eq.  75 
gives  now 


J^(jA  (_a£ — \wl ,  [W_Bg  (1-H)]  zLmu  eL 

z*  Vc  /  Vi-T  n  /w  p  e  p 

f  C  8  r  r 


(81) 


Construction  of  the  response  function  now  proceeds  exactly  as  in 
Section  IV-B,  except  that  Eq.  77  is  used  in  place  of  Eq.  59,  and  Wg 
replaces  W.  The  final  result  is: 


m*/m 

p’/p 


AWg  (1  +  8  *n  Cf)  +  ^  ns  +  ng  (X-l) 

X  +  ^  E  ~  (1+28  tn  I,)  -HA  +  -*■  -  1 

A  C  p  £  C 


which  again  has  the  form  of  Eq.  71,  with  B  now  defined  as 


W  c  n 

B  -  —  (1  +  8  la  Zf)  +  -2-  -5- 
n  f  c  An 


(82) 


5  Note  that  because  fhe  ”1"  has  been  dropped, 
the  solid  is  a  small  fraction  of  the  heat  released 
also  one  is  still  restricted  to  H  >  0  -  of  remarks 


the  heat  transfer  to 
in  the  gas  phase,  and 
following  Eq.  58. 
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t 

I 


The  correct  limit  for  u  *  0  implies  now 

_  c 

AW$  (1  +  8  In  cf)  +  ^  ns  -  n 


A(l-H)  +  ^  E  ^ 


(1  +  20  £n  C.)  +  -E- 

f  c 


(83) 


To  check  that  thi3  is  identically  satisfied,  substitute  Eq.  73  into 
Eq.  11  to  find 


E  -  A  (1-H) 
7  0 
cf 


(84) 


It  is  a  simple  matter  of  substitution  to  verify  that  the  condition  Eq. 
83  is  in  fact  satisfied. 


For  8*0,  the  results  of  Ref.  2  are  of  course  recovered.  Moreover, 
since  the  response  function,  as  a  function  of  frequency,  has  always  the 
same  form,  Eq.  71,  the  same  numerical  results  are  found,  although  with 
different  values  of  some  of  the  parameters  if  different  formulas  B  nre 
used.  A  convenient  way  of  comparing  is  to  define  a  new  quantity  Hg  so 
that  B  here  has  the  same  dependence  on  Rg  as  B  of  Eq.  72  has  on  H,  i.e.. 


;W  (1+  s  & 

n  p  re  An  n 


2n  (1-Hg)  + 


_E 

c 


n-n 


(__8) 


Since  Wg  is  given  by  Eq.  81,  one  finds  the  relation  between  H  and 

1-Hg  -  (1-H)  (1  +  6  in  Cf) 
or 

Hg  -  H  (1  +  B  in  ?f)  -  0  £n  ;f  (85) 

The  meaning  of  Hg  is  thwt  where  H  is  used  in  Ref.  2,  Hg  should  appear 
to  account  for  a  combustion  zone  displaced  from  the  surface.  That  is, 
the  same  numerical  results  are  obtained  here  for  a  value  of  Hg  *  H0, 
say,  as  are  obtained  in  Ref.  2  for  H  -  Hc.  But  if  the  calculations  are 
Interpreted  in  terms  of  Qs,  different  results  are  obtained.  In  Ref.  2, 
Qs  *  c(Ts-Tc)H,  whereas  if  Hg  is  used,  Eq.  85  shows  that 
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Nov  Cf  is  large,  so  that  it  is  quite  possible  that  $  In  is  a 
significant  correction.  Indeed,  if  6  la  »  1,  then  Hp  «  8  in  Cf(H-l) 
which  is  not  only  large,  but  Hg  has  a  different  sign  from  K.  Since  it  is 
clearly  not  precise  to  assume  that  combustion  really  begins  exactly  at 
the  Interface,  it  appears  that  the  modification  suggested  here  may  be 
important  not  only  quantitatively  but  also  qualitatively. 

For  example,  in  Ref.  9,  the  effect  of  aluminum  on  the  transient 
response  has  been  examined  by  relating  changes  of  H  to  the  heat  capacity 
of  molten  aluminum  collecting  on  the  surface.  Thus,  an  increase  of  H, 
due  to  the  heat  absorbed  by  the  aluminum,  accompanies  greater  accumula¬ 
tion  of  the  metal  on  the  surface.  It  was  chen  found  by  numerical  re¬ 
sults  that  increasing  H  (i.e.,  more  aluminum  on  the  surface)  is  a  stabil¬ 
izing  effect:  it  tends  to  reduce  the  peak  in  the  response  function.  This 
was  offered  as  an  explanation  for  the  observed  effect  that  the  addition 
of  aluminum  to  a  propellant  tends  to  suppress  pressure  oscillations  in  a 
chamber . 

On  the  other  hand,  in  both  Ref.  4  and  5,  an  increase  in  H,  or  the 
corresponding  parameter,  was  found  to  enhance  the  peak  in  the  response 
function.6  The  reason  for  the  difference  between  these  conclusions 
seems  to  be  that  account  is  taken,  in  the  last  two  works,  of  the  thermal 
wave  in  the  surface  layer.  This  tends  to  enhance  the  peak  in  the  same 
way  that  the  thermal  wave  in  the  solid  causes  the  peak  originally  for  a 
pure  solid.  If  the  thickness  of  aluminum  is  comparable  to  the  wavelength 
of  the  thermal  wave — a  likely  circumstance — this  influence  cannot  be  ig¬ 
nored.  Hence,  the  conclusions  of  Ref.  9  must,  for  the  present  at  least, 
remain  qualified  in  an  Important  respect. 

A  remark  should  be  made  in  regard  to  the  interpretation  of  the  in¬ 
fluence  of  H  (i.e.,  Qs) -  By  combining  Eq.  21  and  25,  the  overall  energy 
balance  for  the  conversion  of  cold  solid  to  gaseous  reaction  products  is, 
in  the  steady  state. 
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Thus,  when  different  numerical_results  for  the  admittance  function  are 
obtained  by  changing  H  *  Qs/c(T3-Tc),  this  is  necessarily  accompanied 
by  changes  in  the  values  of  Qf,  Tf,  or  Ts  so  that  Eq.  87  is  satisfied. 
This  is  not  troublesome,  but  must  be  recognized  in  the  interpretation 
of  parametric  studies. 

D.  Quasi-Static  Results  by  Integration 

of  the  Linearized  Differential  Equations 

Certain  aspects  of  the  problem  are  clarified  by  finding  the  per¬ 
turbations  by  integrating  the  differential  equations.  It  is  also  neces¬ 
sary  to  follow  this  route  if  one  wishes  to  relax  the  assumption  of  quasi¬ 
static  behavior.  When  the  time  derivatives  are  ignored  in  Eq.  52  and  53, 
and  the  variables  are  written  as  T  ■  T  +  T',  etc.,  one  finds 

3m' /3x  -  0  (88) 


-  3T' 


me 

P 


3x 


Qfw' 


+  m'c 

P 


dT 

dx 


(89) 


In  order  to  simplify  the  computations,  only  the  limiting  case  of  com¬ 
bustion  extending  to  the  solid-gas  interface  will  be  considered  here. 

It  is  also  perhaps  less  confusing  for  the  time  being  to  use  the  position 
variable  x  rather  than  the  dimensionless  £.  The  solutions  for  the 
steady-state  heat  transfer  and  temperature  are 


q(x)  -  kg  ~  =  Qf  w  l  [1  -  e-0]  (90) 

QfW*2  a 

T(x)  -  Tf  -  - -  [6  -  1  +  e  H]  (91) 

8 


where  L  *  kg/mCp  and  0  *  (xf-x)/£  and  the  boundary  condition  qf  ■  0  at 
x  ■  Xf  has  been  used. 

Equation  88  of  course  gives  the  result  that  m*  is  constant  through 
the  gas  phase.  It  is  a  straightforward  matter  to  integrate  Eq.  89  to 
find,  for  the  case  when  w'  is  uniform  in  the  gas  phase: 

<’  "  kg  f1  "  *l'  ^  +  ^  W’£(1-e"6) 

+  Z"  Qf  vt  [(6+l)e"9  -1]  (92) 

m 
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(93) 


T’  -  T-’  -  ~  (l-e“0)  -  Qf  w’  [9-1  +  «  1 

f  f  g  g 

-  zr  V  T  [2(1_e”0)  “  5  <1+e“®)l 

®  g  I 

where  5  ■  mcp(xr-x)/k„,  end  t  has  m  in  place  of  m-  The  integrations  have 
been  performed  from_tne  surface  (x  ■  0)  to  the  mean  position  of  the  outer 
edge  of  the  flame  (xf ) ;  extension  of  the  range  to  x  -  Xf  +  xf ' ,  the 
instantaneous  position,  merely  adds  terms  of  higher  order  in  small  quan¬ 
tities.  Thus,  Tj ’  and  qj'  denote  fluctuations  at  the  mean  positio  ;  at 
the  instantaneous  position,  Tj'  +  TjT  and  qj'  ■  0,  but  note  that  qf'  is 
nonzero.  The  situation  is  sketched  in  Fig.  3. 


FIG.  3.  Conditions  at  the  Downstream  Boundary  of  the  Flame  for 
Quasi-Static  Behavior. 
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Now  if  Eq.  90  and  91  are  linearized,  the  perturbations  oust  be 
exactly  Eq.  92  and  93.  but  two  points  should  be  noted.  First,  there 
are  contributions  from  fluctuations  in  6,  for 


6  - 


e  +  0'  - 

Hence,  in  the 
6 

e  1 


me  me  xr 

(x  -  x)  »  — E-LI. 
k  uf  }  fc 
g  g 


me  x. 

B-J. 


_i  xf 

lr! -  +  JL 


m  x 


f  J 


me  me  x,  / 

r6  wt-x)  +  -H  ( 


exponentials, 

0+0'  0  0' 
e  «  e  e 


e6  (1+0') 


Second,  since  Eq.  90  and  91  extend  all  the  way  to  the  edge  of  the  flame, 
perturbations  at  the  downstream  edge  refer  to  the  instantaneous  and  not 
the  mean  position  of  the  flame.  Thus,  Tj  in  Eq.  91  gives  Tf*,  not  Tj' , 
and  no  term  in  qj ' ,  which  appears  in  Eq.  92  and  93,  is  found.  By  ex¬ 
pansion  in  Taylor  series,  similar  to  Eq.  6  and  7,  one  has 


(94) 


(95) 


(96) 


When  these  relations  are  taken  into  account,  the  desired  agreement  is 
found.  It  may  be  noted  that  since  the  mean  heat  transfer  vanishes  at 
the  flame  edge,  (dT/dx)J  ■»  0  and  to  first  order  Tf'  *  TJ' ,  but  the  heat 
transfer  fluctuation  at  the  mean  position  does  not  vanish  to  first  order. 
The  fluctuation  of  heat  q~’  is  therefore  due  to  the  fluctuation  of 
energy  release  associated  with  the  change  of  flame  thickness. 
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The  last  point  merits  closer  attention.  Consider  the  energy 
balance  for  the  fluctuations  in  the  region  xf  £  x  <_  Xf’  shown  in  Fig.  3. 
One  has  for  Cp  constant 

Qfwxf’  +  Cp  (S rf’  +  o'T£')  -  ep(iTf’  -  m'Tf)  -  qf*  - 
or 

<Tf'  -  <tf’  +  ocp  (Y  -  Tf ')  +  m'cp  (T.  -  Ty> 


+  Qfwxf'  **  Qfwxf  (97) 

since  both  temperature  differences  vanish,  to  fir'*  order,  and  q/  *  0 
by  definition  of  the  outer  edge  of  the  flaue.  It  follows  from  the 
differential  equation  for  the  mean  temperature  t:v;t  kgfd^T/dx^)^  *  -QfW 
and  2q.  96  and  97  are  the  same.  This  exercise  shows  now  it  Is  that 
even  though  conditions  downstream  of  the  flame  are  uniform,  the  f  3  me 
temperature  can  fluctuate,  giving  rise  to  the  entropy  waves  discussed 
in  Ref.  2. 

When  evaluated  at  x  -  0,  Eq.  92  gives  the  expression  for  heat 
transfer  to  the  surface,  identical  to  that  discussed  In  Section  1V-A. 

By  use  of  Eq.  92,  many  of  the  terms  in  Eq.  93  can  be  identified  as 
q'(x)  so  that  Eq.  93  may  be  written 

/  /  _ 

T’  -  V  -«Jf  k +1  k‘w’Qfre 
8  8 

2 

+  rQf  Ue-D  +  e'9)  (98) 

m  g 

Evaluated  at  x  ■  0,  this  gives  the  perturbation  of  the  overall  energy 
balance  for  the  gas  phase: 

v' 5  v'  ■  Qf  <**{' +  «f>  -  "Vh  •  V 


-  5c  (T/  -  T  ')  (99) 
p  f  s 

Another  way  of  arriving  at  the  same  formula  is  to  consider  a  control 
volume  extending  from  x  ■  Q  to  x  «  Xj  +  x,1  or  integrating  the  energy 
equation  directly  over  the  entire  gas  phase: 
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X 

d_ 

me  T 

-  k  ? 

dx  = 

dx 

P 

g  dx 

Qf  v  dx 


With  T  *=  T  +  T',  etc.,  and  retaining  only  terms  of  first  order  in 
fluctuations,  this  gives 


T* 


-  k 


ill 

dx 


Xf+Xf 

- 

+ 

m’c  T 

0 

L  P 

Xf+Xf’ 


*  Qf  (xfw*  +  wxf ’) 

and  since  again  (dT'/dx)  _  —  ,  =0,  Eq.  99  is  recovered, 

X  —  X  ^ 

The  computations  above  serve  to  verify  explicitly  for  the  special 
model  of  uniform  combustion  that  the  procedure  of  Section  IV-A  based  on 
the  steady-state  solution  does  indeed  give  exactly  the  same  results  ob¬ 
tained  by  integrating  the  differential  equations  for  the  perturbations. 

E.  Perturbations  of  Reaction  Rate 
Depending  on  Temperature 

It  was  argued  in  Section  IV-A  that  a  potentially  important  failing 
of  the  model  proposed  by  the  Princ'-'on  group  is  that  the  energy  release 
is  forced  to  be  in  phase  with  the  pressure  everywhere  in  the  gas  phase. 
This  Is  obviously  an  approximation,  but  on  the  other  hand,  treatment  of 
nonuniformities  is  very  difficult,  not  only  because  the  calculations 
rapidly  become  very  involved,  but  also  because  it  is  not  clear  what 
constitutes  a  realistic  representation. 

At  the  present  time,  it  does  not  seem  worthwhile  to  give  up  a  priori 
specification  of  the  combustion  distribution  for  the  steady  state.  And 
in  particular,  if  viewed  in  the  way  suggested  after  Eq.  63,  the  assump¬ 
tion  of  a  uniform  distribution  may  be  accepted  as  a  reasonable  first 
approximation.  However,  this  does  not,  as  argued  previously,  irrevocably 
lead  one  to  use  a  uniform  distribution  of  energy  release  fluctuations 
which  unavoidably  forces  the  energy  release,  or  reaction  rate,  and  pres¬ 
sure  to  be  in  phase. 

The  point  is  that  in  fact  one  has  a  great  deal  of  freedom  in 
respect  to  choice  of  w’;  a  very  simple  alternative  is  offered  here, 
suppose  first  that  the  reaction  rate  is  locally  sensitive  to  temperature, 
such  that  the  steady-state  value  is  approximated  by  a  uniform  distribu¬ 
tion,  but  the  fluctuations  are 
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where  T*  Is  the  local  temperature  fluctuation  and  V  is  a  coefficient 
which  c&n  be  determined  only  by  a  more  careful  examination  of  tne  pro¬ 
cesses  in  the  flame.  It  does  not  appear  possible  to  find  V  as  W  was 
found ,  by  appeal  to  calculations  and  experimental  results  for  the  steady 
state.  This  expression  simply  adds  a  term  to  the  formula  Eq.  81  already 
used.  However,  the  formulas  for  T'  and  q'  in  the  gas  phase  cannot  now 
be  deduced  from  th?  steady-state  solution  which  is  supposed  to  be  fo-  <. 
uniform  distribution  of  combustion.  Rather,  Eq,  100  must  be  substituted 
in  89  and  the  differential  equation  solved.  Incidentally,  it  must 
be  emphasized  that  the  arbitrary  addition  of  the  term  in  Eq.  100  in  no 
way  violates  the  assumption  of  quasi-static  behavior,  so  that  Eq.  88  and 
89  art  still  valid. 


The  computations  required  are  straightforward  but  involved,  and  will 
not  be  covered  here.  A  quick  estimate  of  the  influence  of  this  modifi¬ 
cation  can  re  obtained  in  the  following  way.  Suppose  that  the  influence 
of  temperature  changes  is  most  important  in  the  hotter  regions  of  the 
flame,  and  la  therefore  roughly  represented  by  setting  T*  “  Tj’  in 
Eq.  100.  Calculi *ion  of  the  response  function  now  proceeds  as  in 
Section  TV-C.  Substitution  of  Eq,  100  in  Eq,  77  gives 


(>.8'L  -(■-?)(■• 


_  \  I 


(1  +  26  In  t  )  — 
1  m 


(101) 


__By  use  of  the  perturbation  of  the  ove-all  energy  balance,  Eq.  60, 
Tt'/Ts  can  be  eliminated  from  this  equation,  leading  to 


(102) 


where 
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(1+8  £n  i,^) 


(105) 


The  response  function  is  once  again  found  by  using  the  energy  matching 
condition,  Eq.  18,  the  pyrolysis  law,  Eq.  17,  and  now  Eq.  102  for  the 
heat  transfer  from  the  gas  phase;  the  result  is: 


aVm  ,  ftWBT  (1+6  V  +  Vt)  *  "s1'-1* 

p  '  It  .  c  .  c 

X  +  f  +  -2-  EX-HA  +  -  V  1  -  1 

X  c  \c  T / 

First  checking  the  limit  'u  =  0,  one  has  the  requirement 


(106) 
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f  s  \  c  T ,  1 


c  /  c  , 


which  can  be  rewritten,  after  substitutJon  of  Eq.  103  and  105  as 


AWg  (1+8  Zn  ?f)  ■  n 


Cd  A2 

A(l-H)  +  -R  E  ~  (1+28  In  Cf) 


c  <f6 


+  (n-n  )  D 
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(107) 


39 


*'*$£***-  ’ 


NUC  TP  4668 


where  D  stand*  for  the  denominator  of  Eq.  103-105.  It  follows,  after 
use  of  Eq.  104,  that  D(cp/c  -  V<r)  *  cp/c  and  the  condition,  Eq.  107,  is 
exactly  Eq.  83,  which  has  already  been  found  to  be  satisfied  identically; 
thus,  the  limit  for  a  ■  0  is  correctly  met  for  any  value  of  V. 

Naturally,  the  response  function  again  has  the  fo  m  of  Eq.  71,  but 
with  B  given  by 


n  / 

-2-  1  - 
n.  \ 


£_  v 
c  VT , 
P 


(108) 


As  in  Eq.  82  and  85,  the  results  can  be  interpreted  in  terms  of  a  new 
quantity  Hgj,  defined  so  that  B  here  has  the  same  dependence  on  Hgf  as 
B  had  on  H  in  Eq.  70;  by  use  of  the  definitions  of  Eq.  103,  82,  and  81, 
ia  defined  by  the  equality 


2n(l-H)  + 


“  (1+8  in  C{) 


i?i]- 


82(1-H) 


which  eventually  gives 


1  +  f-  V(1  -s£)  (1+  bin  Cf> 


c  n  . 

f_£_£  !  .I 

c  n  \  c. 


.  i  .  c  ,  n-n  .1  c  n 

vi)4[2nKi)  +  ^(  - 


1+8  in  C,  a  -  8  £n 

V  - H  -mr1  +  —I -+— 1 


0-09) 


where 


V  (l  -  ~ )  <1+&  sf> 


(110) 


Equation  85  is  of  course  recovered  for  V  *  0. 


It  is  not  obvious  how  to  evaluate  V — although  for  a  genuine  diffu¬ 
sion  flame  one  might  try  using  the  temperature  dependence  of  the  diffu¬ 
sion  coefficient — so  to  assess  the  influence  of  this  correction,  suppose 
once  again  that  8  in  Cf  >>  1.  Then  Eq.  109  becomes 
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H-l 


+  1 


(111) 


It  appears  that  this  correction  is  not  quite  so  large  as  that  shown  in 
Eq.  85  for  8  in  Sf  >>  1.  When  V  is  nonzero,  the  influence  of  a  dis¬ 
placed  combustion  zone  is  in  effect  blunted.  That  this  is  the  case  may 
be  interpreted  roughly  as  follows.  If  combustion  begins  downstream  of 
the  interface,  then  the  heat  transfer  to  the  surface  is  less  than  that 
when  combustion  begins  at  the  surface.  On  the  other  hand,  for  V  0, 
the  heat  transfer  is  increased  according  to  Eq.  102.  Hence,  the  two 
effects  tend  to  compensate  one  another. 

Admittedly,  the  preceding  calculations  are  approximate.  The  intent 
is  to  demonstrate  that  the  results  of  Ref.  2  and  subsequent  applications 
in  Ref.  9  and  10  must  be  qualified  in  respects  which  are  not  considered 
in  those  works. 
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In  spite  of  the  assumptions  involved,  the  model  based  on  uniform 
combustion  is  nevertheless  a  most  convenient  one,  and  seems  to  be  an 
acceptable  zeroth  approximation  to  whe  real  problem.  It  is  particularly 
useful  for  studying  an  aspect  of  the  general  p  oblem  whi~h  has  received 
too  little  attention,  namely,  behavior  when  the  assumption  of  strict 
quasi-static  is  not  erforced.  The  only  treatment  of  this  question  is 
Ref.  t,  which  is  not  only  complicated  but  is  based  on  a  very  thin  flame 
front  as  the  model  for  the  gas  phase. 

The  two  conservation  equations,  Eq.  52  and  53,  will  be  solved  with 
the  time  derivative  terms  nonzero.  By  use  of  Eq.  52,  the  left-hand  side 
of  Eq.  53  can  be  written 


c 

P 


fc  <*>  +  cp 


3_ 

3x 


(mT) 


and  with  the  perfect  gas  law,  the  first  term  is  (c_/R)  3p/3t  which  can 
be  grouped  with  3p/3t  on  the  right-hand  side.  It  is  best  to  work  with 
dimensionless  variables;  in  terms  of  £  defined  earlier,  the  equations 
can  be  written 


(11?) 


3C<  +  5  3l  [t(1'w)] 


(113) 


where  t  ■  T/T8,  y  ■  m/m,  and  t  now  stands  for  the  dimensionless  time, 
the  scale  being  pTg  kg/5rcp,  i.e.,  t  stands  for  (m2cp/pgk  )t.  In  order 
to  get  results  without  excessive  complications,  it  will  Be  supposed  that 
w  is  always  independent  of  position.  This  means  that  dependence  of 
reaction  rate  on  temperature  can  be  accounted  for  only  in  the  somewhat 
artificial  manner  discussed  in  Section  IV-D. 


When  the  quasi-static  assumption  is  dropped,  it  is  no  longer  possi¬ 
ble  to  obtain  neat  closed  solutions.  Several  methods  of  approaching  the 
problem  will  be  indicated,  but  eventually  only  the  first-order  correction 
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for  the  case  of  harmonic  oscillations  will  be  treated.  With  w  dependent 
only  on  time,  Eq.  113  can  be  integrated  once--divide  by  C  and  integrate 
from  an  arbitrary  position  to  the  flame  edge  where  in  all  cases 
dr/dc  >»  0.  After  integration  of  the  first  term  by  parts  and  combination 
of  the  various  pieces,  one  finds 


(114) 


Thus,  Eq.  112  and  114  constitute  a  coupled  pair  of  firso-order  equations 
not  yet  restricted  to  small  motions.  For  the  purposes  of  this  section 
only  small  fluctuations  are  considered;  the  linearized  forms  of  Eq.  112, 
113,  and  114  are 


1  3x ' 


W  at" 


(115) 


V  3  -  Id  lp'\  .2  w' 
7T  -  c  TT  ^p’)  -  v  7T  r*  “  A  5“ 


(116) 


(117) 


The  fact  that  TT  =  1  has  been  used,  and  pj'  stands  for  the  fluctuation  of 
mass  flux  at  the  flame  edge.  As  before,  the  fluctuation  of  flame  thick¬ 
ness  in  Eq.  117  must  be  found  from  its  connection,  Eq.  29,  with  w  and  m. 
Equation  117  is  identical  with  Eq.  98;  to_show  this  l.c  is  necessary  to 
replace  qf'  by  use  of  Eq.  96,  and  to  use  8  rather  than  C  as  the  inde¬ 
pendent  variable. 

Second-order  equations  for  both  p'  and  x '  separately  can  be  con¬ 
structed.  For  example,  solve  Eq.  117  for  p'  and  substitute  into  Eq.  115 
to  find 


43 


NWC  TP  4668 


3  [  1  9tM  1  3t’  /C  dr  \ 
3C  \T  3C  /  "  7*  3t  (t7  i C  i 


.(S-  " 

'  r2  dc  / 

* 


+  t.'  +  (A2  -tn 


-1 

5f  J 


(118) 


Solution  of  this  equation  will  give  t*  and  then  Eq.  108  can  be  used  to 
find  y'  by  direct  integration  from  the  surface  outwards: 


rt 


v  + 


c 


1 

5? 


ill  .  i.  1.  (£l)l 

3t  Ct  dt  Vp  I  \ 


dC 


(119) 


Note  that  y8'  is  connected  to  the  surface  temperature  fluctuation  ts*  by 
the  pyrolysis  law,  so  that  Eq.  119  can  be  used  as  a  formula  for  y*',  the 
quantity  ultimately  required.  Since  t  is  an  awkward  function  of 
analytical  solution  to  Eq.  118  does  not  saem  possible,  although  Eq.  118 
and  119  could  be  used  as  the  basis  for  numerical  calculations.  An 
interesting  approximation — at  first  sight  extremely  crude — is  to  set  the 
mean  temperature  equal  to  some  average  value  ta  everywhere  so  that 
dT/dc  -  0  as  well;  then  Eq.  118  is 


.2  »V  .  3V  .  1  d_  (  gll  2 

3C*  T  3t  Y  dt  \  p  /  w 


(120) 


Since  ?2  32t’/3C2  •  (kg/mcp)Z  (racp  3t’/3x  -  kg  32t*/3x2]  (of  Eq.  32), 
the  left-hand  side  is  the  ordinary  heat  conduction  equation  with  con¬ 
vection  of  energy  by  the  mean  flow.  The  right-hand  side  is  of  course 
the  fluctuation  of  energy  release  associated  with  the  chemical  reactions. 


To  see  how  close  Eq.  120  is  to  a  more  acceptable  approximation, 
substitute  Eq.  115  for  C3y'/3C  in  Eq.  116: 


32i 

i? 


1 

3t  ' 

ft Hi) 

fL/£ 

1)  +  A2  + 

U*L\ 

T 

dt 

l  Y  / 

dt  V- 
p 

'  w  - 

r  dc  1 

(121) 


This  is  exact  within  the  initial  assumptions  used,  but  the  difficulty  of 
course  ia  that  y'  on  the  right-hand  side  is  not  known.  In  addition  to 
the  terms  in  Eq.  120,  there  appears  on  the  right-hand  side  a  term 
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d/dt(p'/p)  which  is  a  consequence  of  compressibility,  and  the  term  with 
y  in  it,  which  represents  transport  of  mean  energy  by  fluctuations  of 
the  mass  flux.  A  more  reasonable  approximation  is  to  set  t  ■  7a  on  the 
left-hand  side  of  Eq.  121,  retain  dx/di;  on  the  right,  and  then  worrv 
about  what  to  do  with  p'.  This  procedure  will  be  followed  for  the  case 
of  harmonic  oscillations  whose  frequencies  are  not  too  high. 

A.  Solution  to  the  Gas  Phase  for 
"Low"  Frequency  Oscillations 


The  time  derivatives  are  now  replaced  according  to 


3_ 

at 


ifte 


where  Q  =  <cui/r2  is  the  dimensionless  frequency  arising  in  the  problem 
of  unsteady  heat  conduction  in  the  gas  phase,  and 


e 


(122) 


_2 

is  a  small  quantity  of  order  10  .  It  is  a  little  easier  to  work  with 

the  variable  C  :  C  -  e*  so  that  Eq.  (121)  and  (116)  are 


d2t  ’  drj_  ige  , 
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(123) 


dp 1  _  ifle 

d?  ~  7 


(124) 


In  later  calculations  it  will  be  assumed  that  0e  is  small,  the 
quasi-static  assumption  of  course  corresponding  to  fie  =  0.  To  interpret 
the  meaning  of  this  parameter,  rewrite  it  in  the  form 


fie 


P  k 

■4*-* 
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The  ratio  of  scan  heat  transfer  rate  at  the  surface  to  the  rate  of 
•nergy  transferred  by  convection  appears  in  brackets;  this  is  multiplied 
by  the  ratio  of  the  average  transit  time  (xf/va)  for  an  element  of  gas, 
to  the  period  l/w  of  the  oscillations.  Thus,  it  is  the  product  of  these 
which  matters,  not  the  relative  importance  of  heat  conduction  or  transit 
time  alone. 

The  relatively  weak  assumption  will  be  made  that  7  *  Ta  in  ifle/t; 
in  any  case,  7  varies  smoothly  from  unity  at  the  surface  to  a  value  of 
three  or  four  at  the  flame,  which  is  a  very  small  change  relative  to  c2 
in  the  second  derivative  of  Eq.  121,  i.e.,  compared  with  the  first  two 
terms  in  Eq.  17.  It  appears  that  the  gross  behavior  of  the  results 
should  not  be  greatly  affected  by  this  approximation,  and  by  doing  so 
one  has  the  opportunity  to  obtain  results  quite  easily.  For  now  Eq.  123 
is  formally  integrable.  The  solutions  to  the  homogeneous  equation  are 
exp(X+£)»  exp  (A_£),  where 
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Thus,  some  of  the  integrals  in  the  solution  for  t'  can  be  carried  out 
and  r '  is 


M  x-e 
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(128) 


The  :onstants  and  C2  are  set  by  requiring  that  t'  =  if’  =  tf'  and 

(3t’/3£)—  *  q— '/me  T  **  (Q,w/mc  T  )x,'  (see  Eq.  96  and  97)  one  finds: 
ffps  t  p  s  r 
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' 2  X_-X+ 


-  X+  (Tf'-K)  +  z 


¥ 


me  T 
P  s 


(130) 


With  these  expressions  for  the  constants,  the  solution  for  r’  can  be 
written 
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=*  K  + 
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Since  y'  is  given  by  Eq.  119,  Eq.  131  is  still  not  an  explicit  solution, 
but  it  is  a  convenient  form  for  iterative  approximations,  from  which  both 
x*  and  (dx'/d£)  .  can  be  computed. 

Ot 

The  zeroth  approximation  is  found  by  setting  fis  m  0  in  both  Eq.  119 
and  131;  this  is  of  courre  the  quasl~static  case  considered  in  Section  IV: 
Eq.  119  gives  y'  a  constant  through  the  gas  phase,  and  Eq.  131  reduces 
to  Eq.  93  or  98  if  prope:  account  is  taken  of  the  variables,  0  *  Cf~5* 

The  next  and  higher  approximations  may  be  found  by  expanding  in  powers  of 
ifie.  Only  the  correction  or  order  fie  V7ill  be  considered  here,  so  that 
the  results  are  valid  orly  for  relatively  low  frequencies.  From  the 
definitions,  Eq.  125  and  128,  it  follows  that  X+X_  ■  -ifie/7a  exactly, 
so  that 
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xa 

ifie 


and  to  first  order 
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With  these  used  in  Eq.  131,  one  eventually  finds 
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which  agrees  exactly  with  Eq.  93  since  A2  *  (QfkgW/iPc  ^g) ;  the 
constants  o^,  . ..,  6j,  are  defined  in  the  nomenclature. 

The  first-order  correction  to  the  fluctuation  of  mass  flux  is 
given  by  the  integral  of  Eq.  124: 


„.<!>  .  „  .(1)  + 


p  J 
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The  mean  temperature  profile  as  a  function  of  £  is 


[(?f-«>  -  (1-e  ^  S] 


T  *  Tf  -  A2  5  'r  /f  "  £' 


which  prohibits  easy  evaluation  of  the  integrals.  Hence,  in  the  interests 
again  of  obtaining  results,  7  is  approximated  by  7a.  With  Eq.  126,  one 
eventually  finds  for  y'(D: 


+1  k-T.tjc  +  Bj 


+  yx  {e*  (£-l)-l}  +  f 


(134) 


Note  that  the  mean  value  £s  «  0  has  been  used  since  fluctuations  £s ' 
contribute  terms  of  second  order  only.  Equation  134  is  required  in  the 
first-order  correction  t*'*'  for  the  temperature  fluctuations.  However, 
the  entire  profile  t*'*'  is  not  required  for  the  computation  of  the  re¬ 
sponse  function.  One  needs  only  the  gradient  at  the  surface,  which  from 
Eq.  131  is 
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Henceforth,  only  first-order  terms  in  fie  will  be  considered,  and  <  lu~ 
stion  of  Eq.  135  is  relatively  easy.  First  in  the  integral 

>5f  _  (?f 

g’df  (X+e  '  X-e  )d*»  *  (y'(0>  +  i^w,(1))f|(e"?  +  |S£)d5  => 


i * v°) a2  l  -  (^  +  l)e  +  ine  Ks - (t,  -  x  )  +  p,(1)  —e~^  d? 


—  V  I  c  v  / 

U  f  8 


With  the  very  good  approximation  that  terms  containing  exp(-C^)  can  be 
neglected,  one  finds 
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If  the.  remaining  terms  are  expanded  to  order  ifle  and  terns  of  order 
expC-iff)  are  dropped.  Eq.  128  becomes 
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The  zeroth  order  terms  give 
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(136) 


which  la  exactly  the  qua  -static  results,  Eq.  92,  evaluated  at  the 
surface.  The  first-order  .erms  give 
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The  definitions  of  sj.  and  6].  can  be  substituted  to  give  a  more  explicit 
form.  In  doing  s Oj_  further  terms  are  neglected,  based  on  the  fact  that 
A2  *  (Qf/c«¥8)(£n  Cf)~l  *  Qf /(CpTs5f )_  is  much  larger  than  unity  (cf. 
remarks  following  Eq.  57.  Although  Zf  »  mCnXj/k-  may  be  10  -  20, 
Qj/cpTg  may  be  as  large  as  5  -  6,  so  that  the  assumption  A2  «  1  may 


-! 


NWC  TP  4668 


sometimes  be  borderline — it  should  be  checked  in  individual  cases* 
However,  it  i6  true  that  the  terms  containing  A2  will  make  a  smaller 
contribution  in  any  case,  and  if  they  are  retained,  subsequent  manipu¬ 
lations  become  very  tedious  indeed.  It  appears  that  the  major  behavior 
will  be  accounted  for,  and  hence  the  essential  implications  of  the  cal¬ 
culation  will  be  found,  if  the  approximation  is  made.  The  relative 
simplicity  of  the  final  formula  F.q.  165  partly  justifies  this  step. 

Then  Eq.  137  reduces  to 
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Some  terms  containing  A2  have  been  retained  to  show  the  only  first- 
order  contributions  on  the  right-hand  side.  However,  these  may  be 
neglected  compared  with  the  remaining  terms,  and  one  is  left  with  the 
quite  simple  result 
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There  are  no  first-order  fluctuations  appearing  on  the  right-hand  side, 
which  makes  later  calculations  very  much  shorter  than  if  Eq.  138  is  used. 
This  i^s  of  course  a  consequence  of  the  fact  that  the  normalized  thick¬ 
ness  Cf  is  very  much  greater  than  unity.  It  is  by  no  means  obvious  that 
all  of  the  labor  following  Eq.  131  would  culminate  in  such  a  result  as 
Eq.  139. 

A  second  necessary  relation  among  the  heat  transfer,  surface  temper¬ 
ature,  and  flame  temperature  fluctuations  can  be  deduced  by  direct  inte¬ 
gration  of  the  linearized  energy  equation,  Eq.  116^_  After  change  cf 
variable  from  Q  to  5  and  integrating  over  0  <_  £  £  Cf ,  one  finds 
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Once  again  expand  all  quantities  to  first  order  in  fle  to  find  the  zeroth 
and  first  order  pieces  of  Eq.  140: 
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Note  that  because  of  the  factor  Kf,  the  last  term  is  dominant.  Equation 
141  is  of  course  the  same  as  Eq.  99;  to  show  this  explicitly  requires 
use  of  Eq.  30  which,  when  linearized  and  expanded,  gives 


w^o)  ,  y.(o)  _  V 


(143) 


for  the  quasi~static  contributions.  Recall  also  from  Eq.  97  (which  is 
valid  for  all  linear  motions)  that 


P  « 

p  •  r 


(144a) 


(144b) 


Equations  141  and  133  may  be  used  to  find  the  fluctuations  of  flame 
temperature,  xj* . 
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The  only  task  remaining  so  far  as  the  gas  phase  is  concerned  is  to 
find  formulas  for  the  fluctuations  of  the  flame  thickness.  For  the 
quasi-static  problem,  Eq.  30,  w  *  m/xf,  can  be  used;  thus  giving  Eq.  142 
but  in  the  case  of  nonquasi-static  motions,  one  tuat  return  to  the  origin 
of  the  term  Qjw.  If  one  examines  the  balance  of  enthalpy  for  the  ga6 
phase,  by  use  of  a  control  volume,  one  finds  by  comparison  with  the  inte¬ 
gral  of  the  unsteady  energy  equation  that  the  following  equality  must  be 
met: 
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where  hf_,  hsr  are  reference  enthalpies.  This  corresponds  to  Eq.  28  for 
steady  flow;  but  here,  instead  of  Eq.  30,  one  has 


w  * 
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Hence,  th»  perturbations  are  related  by 


or,  for  the  first-order  terms: 


(143) 


The  value  of  hfr  depends  on  the  thermodynamics  of  the  reactions  in  the 
gas  phase.  If  the  gases  leaving  the  surface  burn  completely  to  form  the 
elementary  products  in  terms  of  which  the  reference  enthalpy  is  defined, 
then  hfr  *  0  and  Qf  *  hgr,  Rather  than  complicate  the  results  with 
another  parameter,  this  limiting  case  will  be  supposed  here  so  that 


(146) 
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which  five*  Cf'^  '  tor  use  In  Eq.  144b.  Because  of  ths  fora  of  Eq.  139, 
this  result  Is  not  required  for  celculetion  of  the  response  function; 
however,  it  iHist  be  used  in  Eq.  142  to  find  the  fluctuations  of  flame 
temperature.  Kote  that  if  one  assumes  w  to  be  given  by  the  simplest 
form  (Eq.  64),  then  w*  0-)  vanishes,  and  *  Us*^ 

B.  First-Order  Interfacial  Conditions 
* 

The  linearized  forms  of  the  pyrolysis  law,  Eq.  17,  and  conservation 
of  energy  at  the  interface  are  easily  expanded  in  powers  of  fie  to  give 
the  following  relations: 


approximation : 
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Also,  the  linearized  boundary  conditions  at  the  surface  must  be 
expanded  in  powers  13c.  Just  as  in  the  analysis  of  the  solid  phase 
outlined  in  Section  11,  the  solutions  to  the  set3  of  eqwsiions  above  in 
Section  V-A  will  yield  the  values  at  the  origin  x  *  5  -  0,  but  one 
requires  the  values  at  the  Interface  located  at  x(;  again,  for  small 
motions  one  has  Eq.  6  and  7,  and  the  corresponding  formula  for  ys': 
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However,  since  for  the  one-dimensional  problem  <fiT/dx  ■  0,  Eg.  151  shows 
that  yg*  «  u0 *  and  no  correction  is  required  for  surface  motion.  In 
order  to  obtain  the  correct  quasi-static  limit  for  to  ■  0,  it  is  neces¬ 
sary  to  adopt  a  coordinate  system  such  that  the  solid  moves  at  uniform 
speed  "r  plus  a  correction  uQ;  then  conservation  of  mass  applied  to  the 
interface  gives,  instead  of  Eq.  9: 
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In  order  to  have  the  origin  and  the  Interface  coincide  at  all  times  in 
the  limit  w  •  0,  it  is  necessary  that  xg(°)  vanish,  as  argued  in  Section 
III.  Hence,  the  correct  choice  for  u0  is 


(154) 


which  has  the  reasonable  interpretation  that  for  a  small  D.C.  (w  *  0) 
change  of  pressure,  the  solid  is  taken  to  move  from  the  left  at  a  speed 
equal  to  the  original  linear  burning  rate  plus  the  change  of  linear 
burning  rate  corresponding  to  the  change  of  pressure-  Incidentally,  this 
value  for  u0  appears  nowhere  explicitly;  its  presence  here  merely  serves 
as  a  formal  justification  for  the  assertion  in  Section  III  that  values 
at  the  interface  coincide  with  those  at  the  origin  in  the  quasi-static 
limit— T’^  ■  le,  etc-  It  therefore  follows  by  expansion  of  Eq.  6  and  7, 
and  by  use  of  Eq.  153,  that  the  correct  boundary  values  for  the  various 
approximations  are: 
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and  similar  equations  for  the  higher  order  terms. 

C.  First-Order  Correction  to  the  Response  Function 

By  definition,  the  response  function  R  is  the  ratio  evaluated  at 
the  downstream  edge  of  the  flame 
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yf' 

- . .  m - - 

p'/p  p'/p 


R(o)  +  iOe  R(1) 


when  carried  out  tc  first  order  in  iQe.  The  calculations  of  Section  IV 
give  r(°) ,  the  response  function  subject  to  the  quasi-static  approxi¬ 
mation,  R(°)  *  p*  (°)/(p'/p) ;  the  task  here  is  to  compute 


,C1)  ..  _L 


(157) 


Equation  134,  evaluated  at  the  edge  of  the  flame  and  with  terms  of  order 
exp(-Cf)  and  1/Tf  dropped  where  appropriate,  leads  to  the  formula  for 
Pf ( 1> : 
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Thus,  with  Eq.  142,  the  first-order  correction  to  the  response  function 
is 
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Equations  139,  150,  153,  and  156a-b  combine  to  give 
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Now  from  Eq.  64,  with 


Thus,  because  of  the  factor  5f2»  the  dominant  terms  in  Eq.  162  are  the 
terms  involving  w'(°);  approximately,  then: 
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and  with  all  these  simplifications,  the  response  function  to  first  order 
in  fte  is 
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(164) 


It  in  somewhat  surprising  that  such  a  simple  result  can  finally 
be  obtained .  The  primary  rear-n  for  this,  a£  remarked  above,  is  formally 
that  the  normalized  thickness  of  the  flame,  is  large,  and  it  is 
possible  to  throw  away  many  contributions.  The  physical  meaning  of  this 
dimensionless  quantity  may  be  seen  by  writing  it  as 


me 
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The  numerator  is  roughly  the  net  flow  of  enthalpy  through  the  gas  phase, 
and  the  denominator  is  an  estimate^  of  the  heat  transfer  from  the  gas 
phase  to  the  svlid.  Thus,  large  £f  implies  a  relatively  small  heat 
transfer  loss  from  the  gas  phase  to  the  solid,  a  connection  which  also 
follows  from  Eq.  49  and  other  results  arising  in  the  solution  to  the 
ateady-stace  problem. 


Now  the  second  term  in  the  square  brackets  of  Eq.  164  is,  under 
practical  conditions,  not  huge,  whereas  the  first  term  is;  hence, 
finally,  the  response  function  can  be  written  approximately  as 
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(165) 


The  first-order  correction  in  frequency,  therefore,  mainly  affects  the 
imaginary  part  of  the  response  function. 


The  term  proportional  to  Sf2  in  Eq.  165  arises  from  the  formula  for 
Pf*^),  just  before  Eq.  158.  It  is  effectively  the  integrated  effect 
of  the  energy  release  fluctuation  over  the  gas  phase:  note  that 
w'Otf/w-Tfp’/p  according  to  Eq.  71  and  the  manipulations  following 
Eq.  162.  The  energy  release  enters  through  its  effect  on  the  temper¬ 
ature  profile  in  the  conservation  of  mass,  for  example,  as  shown  in 
Eq.  124  and  subsequent  integrations. 


One  therefore  has  the  following  interpretation  lot  the  major 
correction  shown  in  Eq.  165.  According  to  the  conservation  of  mass, 

Eq.  52,  the  fluctuation  o*  mass  flux  depends  on  the  rate  of  change  of 
temperature  with  time,  and  hence  m'  (i.e.,  y')  is  out  cf  phase  with  T’. 
When  the  entire  gas  phase  is  considered,  the  dominant  contribution  to 
the  fluctuation  of  temperature  is  the  fluctuation  of  energy  release, 
which  is,  in  the  simplest  case,  in  phase  with  the  pressure  fluctuations. 
Bence,  since  the  energy  release  and  associated  temperature  fluctuations 
are  in  phase,  the  net  result  is  a  contribution  to  the  mass  flux  which  is 
out  of  phase  with  the  pressure  changes,  the  last  term  of  Eq.  165. 


The  fact  that  the  correction  term  in  Eq.  165  can  be  quite  large, 
and  yet  the  real  part  of  the  response  function  is  not  much  affected,  is 
a  significant  result.  Earlier  numerical  calculations  in  Ref.  11  indi¬ 
cated  that  the  real  part  of  the  response  function  given  by  the  strictly 
quasi-static  analysis  did  not  differ  much  from  that  given  by  the  results 
of  Ref.  6  which  did  account  for  nonquasi-static  behavior  of  the  gas 
phase.  Moreover,  the  real  part  of  the  quasi-static  formula,  the  first 
ratio  in  Eq.  165,  does  give  remarkably  good  (at  least  qualitatively) 
agreement  with  experimental  results.  However,  the  analysis  of  Ref.  7, 
in  which  both  the  real  and  imaginary  parts  of  the  quasi-static  formula 
were  used  to  interpret  data,  showed  very  serious  discrepancies.  These 
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earlier  observations  are  all  accomodated  qualitatively  by  the  result 
shown  in  Eq.  165,  that  nonquasi-static  behavior  contributes  overwhelm- 
ingly  to  the  "out  of  phase"  component  of  the  response  function. 

It  must  be  emphasized,  however,  that  this  does  not  mean  that  one  is 
free  to  determine  numerical  values  of  A  and  B  by  comparing  only  the  real 
part  of  the  response  with  experimental  results.  In  this  respect,  the 
discussion  of  Ref.  7  is  of  course  valid.  The  charts  developed  in  that 
work  will,  however,  be  substantially  changed  if  the  correction  term  of 
Eq.  165  is  taken  into  account.  This  will  be  treated  in  a  subsequent 
publication. 
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